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ABSTRACT 


There has been considerable speculation as to the energy 
balance of the thermosphere and in particular about the fraction of 
the total energy input supplied by ultra-violet radiation. The 
problem is considerably simplified by considering the energy 
balance of the nighttime hemisphere alone. Sunrise and sxmset 
vapor trail measurements provide data on the wind systems at the 
terminator boundary, and temperature measurements provide 
information on the vertical energy conduction. North- south winds 
from -high latitude vapor trail measurements provide a measure of 
the energy input from auroral processes. The discrepancy between 
the inputs from the day hemisphere and the observed loss rates is 
discussed in terms of ion- neutral processes and gravity wave 
inputs , 
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CHAPTER I 


INTRODUCTION 

1. 1 General Stafcement of the Problem 

The thermosphere has received a considerable amount of 
attention since the advent of satellites and it is therefore 
surprising that at this time there is still considerable uncertainty 
about the relative importance of the various energy inputs and 
transport mechanisms in controlling the energy balance. 

Several studies of the overall energy balance of the 
thermosphere have been performed by Swartz et al. (1972), 
Chandra (1973) and Roble and Dickenson (1973). It may be said 
that these investigations have not yet provided any concrete 
results. While there are good measurements of densities at 
satellite altitudes, there are not nearly as many measurements 
of the temperatures and very few are available in the critical 
lower thermosphere where much of the energy is stored and 
dissipated. Very few rocket vapor trail wind observations are 
available at times other than sunrise and sunset and incoherent 
scatter sounding provides little information on the critical zonal 
winds. While extensive series of measurements have been made 
of the EUV fluxes, there are problems with the absolute calibra- 
tion of the intensities in some of the older measurements obtained 
at higher solar activity levels. The heating efficiency in the 
thermosphere has major uncertainties as well. The Joule heating 
at high latitudes and the processes by which the energy of Joule 
heating is transported to lower latitudes are not well understood. 
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In. addition, the energy transported upwards by tidal motions and 
gravity waves into the thermosphere from below is largely 
■undetermined. 

When the entire thermosphere is considered, the ■uncertainty 
in these sources and sinks are of the same order of magnitude 
as the EUV heating which is generally believed to be the major 
heat input and it is difficult to determine whether another major 
source is required. It therefore seemed of interest to examine the 
energy balance of the nighttime thermosphere by itself. The EUV 
heating is a very minor input while the loss rates are almost the 
same as those present during the day and so this approach 
provides a good opportunity to study the effects of transport from 
the day hemisphere. In addition, the vapor trail wind measure- 
ments made at sunrise and sunset provide data on the winds at the 
terminator and additional vapor trail measurements at high 
latitudes provide data on the meridional winds at this boiindary. 

It is the general aim of the present study to attempt to 
determine the rate of energy loss from the night hemisphere and 
compare this with the energy inputs from the day hemisphere to 
determine if another major energy source is required. 

1. 2 Previous Related Studies 

1.2.1 Estimates of Solar EUV Radiation 

The importance of accurate estimates of the solar extreme 
ultraviolet (EUV) flux has been recognized for a number of years. 
Decimeter flux measurements (i.e, the ^ index) serve as good 

indicators of overall solar activity, but direct observations of the 
EUV flux responsible for photodissociation and ionization processes 
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in the thermosphere are required for the bulk of all aeronomical 
calculations in this region. The term EUV has been used to 
designate the portion of the spectrum from roughly 1750 ^ to 170 
while XUV, often taken to be equivalent to EUV, includes in 
addition the X-ray spectrum below 170 ^ with no well defined short 
wavelength limit. 

Early low resolution measurements of solar EUV radiation 
were performed by Detwiler et al. (1961) using photographic- 
photometric equipment aboard an Aerobee-high rocket. The results 
were obtained in the altitude range from 148 to 220 km and covered 
the wavelength range from 850 & to 2600 &. A number of factors 
contributed to uncertainty in the results, including a lack of good 
calibration standards and an uncertainty as to the amount of 
atmospheric absorption along the optical path. 

A synopsis of observations of the thermosphere during 1963, 
including estimates of photoionization cross sections, rate 
coefficients and neutral density distributions was presented by 
Hinteregger and Hall (1965). This paper has been most often 
referred to for its compilation of EUV observations by spectro- 
photometric methods in the wavelength range from 1750 X to 1 & . 
Data obtained between July and August of 1963 agreed quite well 
and the iincertainty in these measurements was estimated to be 
less than 30 percent. A much greater uncertainty existed in the 
350 & - 56 X results however, and the fluxes here are thought to 
be uncertain by as much as a factor of 3 to 5. The uncertainty in 
Hinteregger' s results were again largely due to a lack of reliable 
absolute calibration standards in the ultraviolet. 
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A major concern of the EUV measurements has been 
improved estimates of the flux in the Schumann-Runge continuum, 
in the wavelength range from about 1250 ^ to greater than 2000 X . 
These fluxes are responsible for the rapid dissociation of in the 
lower thermosphere and as such, their knowledge has important 
implications for the global distribution of atomic oxygen in this 
height range as well as the EUV energy input. 

Intensity measurements of the solar continuum from 1400 & 
to 1875 & using a rocket borne photoelectric spectrometer were 
made in September of 1968 by Parkinson and Reeves (1969) in an 
effort to improve on spectral resolution and the intensity calibration 
employed in obtaining previous data. The data collected were 
limited to the center of the solar disk. Their results suggested a 
blackbody temperature of approximately 4600 K and the intensity 
measurements at the short wavelength end ('^1400 ^ ) agreed closely 
with the results of Goldberg et al. (1968). The derived intensities 
were lower by a factor of three, however, than the then current 
estimate of Whiting and Purcell (1969). 

The next comprehensive tabulation of solar EUV data was 

published by Hinteregger (1970) and included the region of the 

spectrum from about 1800 ^ to 30 X, The fluxes arrived at here 

-22 -2 

represented medium solar activity (Fj^q y = 130-170 X 10 W m 
Hz and were considerably smaller than those presented in his 
earlier work (Hinteregger and Hall, 1965) and nearly a factor of 
three less than the results of Detwiler (1961) and Parkinson and 
Reeves (1969). This new information, then, suggested a much 
lower intensity of solar flux in the Schumann-Runge continuum 
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with the consequence that previous estimates of dissociation 
rates and EUV heating should be reduced substantially. 

Measurements in the 1300 ^ - 280 S. interval were obtained 
using a scanning spectrophotometer aboard the NASA OSO-III 
satellite while fluxes in the range 280 ^ - 30 X were obtained from 
an AFCRL rocket using an improved scanning monochromator. 

Hinteregger warned against comparing these values with 
earlier estimates for the purpose of deducing solar cycle variations, 
the reason being the poor quality of these earlier estimates and the 
lack of secure absolute calibration in the earlier flights. He 
suggested, rather, that these new estimates supersede the older 
ones and that any knowledge of the variation over a solar cycle 
remains to be determined. 

Roble and Dickinson (1973) calculated global mean temperature 
profiles using as heat sources Hinteregger's (1970) EUV fluxes for 
wavelengths less than 1300 ^ and the results of Widing et al. (1970) 
for wavelengths greater than 1300 Using a heating efficiency of 
33 percent, they balanced these EUV sources against calculated 
values of conduction and 63 micron cooling by atomic oxygen. The 
temperature profiles arrived at were considerably lower than those 
given by the Jacchia, 1971 empirical model. The results suggested 
that the Hinteregger, 1970 fluxes were low and it was found that 
uniformly doubling or tripling these fliixes brings both the calculated 
and neutral electron temperature profiles into good agreement with 
observations . 

Further evidence pointing to a combined uncertainty ih 
thermospheric densities and the Hinteregger, 1970 fluxes was 
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provided by Swartz and Nisbet (1973). This study involved a balance 
between electron recombination rates determined using Arecibo 
electron density profiles and production rates based on Hinteregger' s 
1970 fltixes, both characteristic of medium solar activity. Photo- 
ionization and photo- absorption cross sections were taken from 
Schoen (1969) and Huffman (1969). The results showed that the 
calculated recombination rate exceeds the photoionization rate and 
indicated that an increase by a factor of two or three in the EUV 
flux is required to remove the discrepancy. A serious error in the 
reaction rates was thought unlikely since this would involve all of 
the rates of interest. The conclusion reached was that additional 
energy input to the F-region is required, either by EUV or some 
new source of energy not yet considered. 

In constructing a photochemical model of the daytime F— region 
G. Scialom (1974) used neutral profiles from incoherent scatter and 
rocket borne mass spectrometer results obtained during 1973. As 
an energy source he employed Hinteregger' s 1970 EUV fluxes using 
a scaling factor as a free parameter. Once again, best agreement 
between experimental ionic concentration profiles and model results 
was obtained when the fluxes were increased by a factor of three. 

From a comparison of chemical production and loss rates, 
Prasad and Furman (1974) reached a .somewhat different conclusion 
and questioned the need for increasing observed EUV fluxes. Their 
deduced loss rate profiles were considerably smaller than those of 
the above mentioned authors and in fairly good agreement with 
production rates based on Hinteregger's 1970 fluxes. 

An alternative treatment by Prasad and Furman which used 
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electron temperatures consistent with observed ion profiles 
suggested also that Hinteregger's 1970 fluxes were low by no more 
than 20 percent. To solve the problem, they suggested, as did 
Hinteregger, repeated, reliable measurements of the EUV flux 
under a variety of conditions and/ or an alternative way of 
examining the problem. 

1.2.2 Mid Latitude Wind Observations 

Studies of the thermospheric neutral wind system at low and 
mid latitudes have proceeded along two fronts; the development of 
models of the neutral wind based on measurements of thermospheric 
pressure and temperature distributions, and experimental techniques 
for the direct measurement of wind velocities. 

Using simplifying assumptions appropriate for the thermo- 
sphere, model calculations of thermospheric winds can be 
performed using the equation of motion for the neutral gas 
(Rishbeth, 1972) given by 


dU 

dt 


+ 


2n X u 


F - v^.(U- V) + 


H v^u + 

p 


g 


( 1 . 1 ) 


where U is the velocity of the neutral wind 

2 X U is the Coriolis force ( being the earth's angular 
velocity) 

F is the driving force per unit mass due to pressure 
gradients 

“V ( U - V ) is the so-called ion drag term (v . is the ion- 
ni ® ' ni 

neutral collision frequency and V is the ion velocity. 
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controlled both, by electric and magnetic fields and by 

coupling with the neutral gas) 

2 -* 

(h-/p)^ U is the viscous drag term (p-Zp being the kinematic 
viscosity) 

g is the acceleration of gravity 

Evaluating the ion drag term in the above expression requires 
an estimate of the ion- neutral collision frequency v^. , This 
quantity has been studied by Dalgarno (1964) and Stubbe (1968) and 
has been found to be density dependent and somewhat temperature 
dependent. 

Empirical models of thermospheric pressure distributions 
have been determined (Jacchia, 1965, 1971) based on the results of 
satellite drag measurements. These models display the density and 
pressure distributions as a function of latitude and local solar time 
and show a characteristic enhancement of the pressure, known as 
the "diurnal bulge", near 1400 hours local time. These results 
allow the determination of the force due to pressure gradients 
given by F in the expression above. The problem with this approach 
lies with the pressure gradients in the neutral models. It has been 
shown by Blamont et al. (1973) and Thuillier et al. (1976), that the 
maximum temperature is not in the equatorial region but at high 
latitudes in the summer hemisphere. This difference has a 
profound effect on the meridional pressure gradient as shown by 
Nisbet et al. (1977). 

It has been customary to divide the experimental methods of 
determining upper atmosphere motions into two broad classes 
(Kent, 1970). The first involves the measurement of ion motions, 
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from which the neutral motions are derived. The second involves 
the direct determination of the motion of the neutral gas. 

Ion motions at thermospheric heights are commonly derived 
using the incoherent scatter technique. Here, reflections occur 
from a beam of VHF or UHF radio waves as a result of Thompson 
scattering by electrons and give information on electron concen- 
trations at heights up to greater than 1000 km. The difference in 
frequency between the peak of the thermally broadened' return signal 
and the incident signal gives an indication of the velocity of the 
bulk motion of the ions. 

Using incoherent scatter measurements of electron and ion 
temperatures, and electron density and vertical plasma drift 
obtained at Millstone Hill (42.6 N, 71.5 W) , Roble et al, (1974) 
have calculated neutral winds for two geomagnetically quiet days in 
March and September- of 1970. The derived zonal winds near 300 
km show a general west to east flow during the afternoon, peaking 
near sunset. A similar peak was found to occur in the east to west 
wind near sunrise. The results obtained were based on an electric 
field model derived from quiet condition electric field measure- - 
ments at Millstone Hill (Salah, 1972). As such, the derived winds 
were considered representative of only geomagnetically quiet 
conditions , 

Of the available techniques for direct measurement of neutral 
thermospheric winds, the rocket vapor trail technique has been the 
most widely employed (Kent, 1970; Haerendal et al.,1969). 

The release of several kilograms of sodium at altitudes up to 
200 km has proven an effective means of measuring neutral wind 
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motions near the morning and evening terminators. The sodi\im 
charge is expelled from the rocket and produces a long visible 
yellow trail as a result of resonant scattering of sunlight at 5890^ 
and 5896^, the wavelengths of the sodium doublet. As a result of 
diffusion, trail lifetimes vary from several tens of minutes at 100 
km to roughly 30 seconds at 160 km. An improvement upon the 
sodium vapor technique involves the use of chemiluminescent trails 
which react with atomic oxygen to form luminous trails visible 
throughout the night. Frequently used chemicals are trimethyl- 
aluminum (called TMA) and nitric oxide (Rosenberg and Edwards, 
1964; Rosenberg, 1966). 

Manring et al. (1964) performed an analysis of 22 vapor trail 
measurements performed at Wallops Island, Virginia (37.8 °N, 

75.5 °W), utilizing the sodium vapor technique in the height range 
from 85 to 135 km. Data was analyzed in terms of 24, 12 and 8 
hour components and although large variations of the wind profile 
were observed from flight to flight, a general pattern of wind flow 
was revealed. 

Rosenberg and Edwards (1964) used chemiluminescent trails 
in the altitude range from 90 to 150 km to measure winds 
throughout the night of December 3, 1962 and May 17- 18, 1963 at 
Elgin Air Force Base, Florida (30°N, 87*^W). The results showed a 
consistent clockwise rotation of the wind with height as viewed 
from above. The direction of the wind displayed a sinusoidal 

>a 

pattern with height, the wavelength increasing from 10 km at 100 
km altitude. 

From an analysis of 25 sodium cloud experiments in the 70 
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to 190 km height range and near the morning and evening 
terminators, Kochanski (1964) attempted to resolve the observed 
motions into contributions from (i) the prevailing seasonal drift, 

(ii) tidal components and (iii) internal gravity waves. He showed 
that the winds in general exhibit large variability in both magnitude 
and direction up to about 130 km and a common feature is a rotation 
of the wind direction with height, in agreement with that observed 
by Rosenberg and Edwards (1964). All profiles showed velocity 
oscillations with wavelengths of about 6 km below 120 km altitude 
and increasing to greater than 80 km at 180 km altitude, indicating 
the presence of gravity waves; The analysis suggested that between 
120 and 140 km, during both summer and winter, most of the winds 
near both terminators are in the same direction, indicating that 
either the 12 hour tidal harmonic is dominant, or that tidal 
components are smaller in magnitude than the general drift. 

In similar fashion, Woodrum and Justus (1968) used wind data 
from 38 chemical releases at Elgin AFB (30.4°N, 86.7°W), 
obtained during the period from October, 1962 to November, 1965. 
Their results indicated a predominant vertical wavelength of 23 km 
between 90 and 120 km in agreement with tidal theory (Hines, 1968), 
and a mean downward phase progression from dusk to dawn. Once 
again, the results showed a clockwise rotation of both prevailing 
and diurnal wind components. 

Observations of changes in the orbital inclination of satellites 
allows a coarse determination of the neutral zonal wind velocity in 
a reference frame fixed with respect to the sun. Using this 
principle. King - Hele (1964, 1966) has monitored the orbits of a 
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nximber of satellites and observed changes different from what 
would be observed if the atmosphere rotated strictly at the same 
rate as the earth. Among the results obtained by King-Hele was an 
indication that the atmosphere superrotates in the altitude range 
from 200 to 300 km with an angular velocity of about 1.4 times 
that of the earth. The possible effects of gravitational forces, both 
terrestrial and l\mi- solar were investigated and found to be 
negligible for the purposes of the study. 

Some success has been achieved in deriving wind motions 
from atomic oxygen airglow emissions at both 577.7 nm and 630.0 
nm (Kent, 1970; Nagy et al. , 1971). Successful measurements of 
strong (> 300 m sec meridional winds were obtained at the 
Michigan Airglow Observatory (42°N, 83 °W) during the night of 
March 14 - 15, 1969. In addition, a comparison between airglow 
winds and those deduced from incoherent scatter measurements at 
Millstone Hill during October of 1969 indicated good agreement. 
Important questions remain, however, about the interpretation of 
these observations in deducing neutral gas motions. 

1.2,3 Joule Heating and the High Latitude Neutral Winds 

Joule heating, or resistive heating of the neutral atmosphere 
through collisions between neutrals and ionospheric currents, 
constitutes a large energy source at thermospheric heights, 
particularly at high latitudes. Discussions of the basic parameters 
involved in the theory of electrodynamic heating by the ionospheric 
plasma, including the effects of viscous dissipation, are given by 
Cole (1971) and Egelend (1975). The driving force for the ions 
originates in the ionospheric convection electric field, of magneto- 
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spheric origin. The energy of the field is transferred to the 
neutral atmosphere via collisions between the ions and neutrals. 

This process results in enhanced neutral temperatures and 
pressure gradients which transfer energy from the region in the 
form of neutral winds. 

In order to establish improved global thermospheric models, 
it is essential to understand the effect of high latitude heat sources, 
particularly Joule heating, on the global wind and temperature 
structure of the thermosphere. To this end, a number of authors 
have examined the temporal and spatial variation of the high 
latitude electric and geomagnetic fields and the motions of both 
ions and neutrals at ionospheric heights. 

Neutral wind and velocity measurements in the range from 
90 to 230 km have been made over Esrange, Kiruna (Rees, 1971) 
using rocket borne chemical releases. The results showed that 
high neutral wind velocities in this region result mainly from ion 
drag acceleration. A comparison of disturbed and quiet condition 
results also showed a correlation of the zonal neutral wind 
component above 120 km and the south- north magnetic perturbation 
vector. This relationship arises since ionospheric currents and 
neutral winds are both functions of the electric field and the 
electron density. Under disturbed conditions, high neutral wind 
speeds have been found near 130 - 140 km and Rees suggested that 
the southerly propagation of these winds gives rise to anomalously 
high winds at lower latitudes following geomagnetic disturbances. 

The relationship between such disturbances and the neutral 
temperature was discussed by Chanin et al. (1972) using 
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temperature measurements from two series of rocket launchings 
during 1969 and 1970. The experiments were performed at Heiss 
Island (80.4°N geographic latitude; 74°N geomagnetic latitude) and 
temperatures were obtained in the range from 120 to 160 km by- 
measuring atmospheric absorption of sunlight through sodium clouds. 
Combined with the results of earlier experiments at this location 
(Polaskov et al. , 1969) the complete data set represented the 
results of approximatel-y 20 launchings, and suggested a relation- 
ship between the neutral temperature and the time rate of change 
of the horizontal magnetic field component several hours prior to 
each measurement. In addition, the data suggested that a slight 
local cooling preceeds the heating from the magnetic disturbance 
in each case. The authors conjectured that most of the observed 
heatings are small scale phenomena resulting from the local 
precipitation of soft (< 1 keV) electrons. The rapid disappearance 
of the sodium clouds during magneticall-y active periods was 
thought to imply a “connection between this localized heating and 
strong winds. 

Using a simplified analysis, Fedder and Banks (1972) 
demonstrated the importance of ion drag at high latitudes in 
establishing the thermospheric wind system. The model employed a 
spatially xiniform time- dependent electric field in the dawn to dusk 
direction, chosen to be representative of the onset of a magnetic 
disturbance, and a uniform magnetic field directed vertically 
downward. The results of this study showed that through ion drag, 
electric fields largely control the flow of neutrals and generate a 
twin- cell pattern over the polar regions. This pattern describes a 
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general night to day flow in the auroral zone and large day to 
night flows over the cap itself where solar tidal winds are the 
dominant driving force. These processes have important implications 
for a global heat balance, since the system is a source of both 
kinetic energy of motion and thermal energy as a result of the gas 
being heated from various polar processes. 

In an attempt to establish statistically appropriate boundaries 
for the polar cap and the auroral zone, Heppner (1972) has analyzed 
the results of OGO-6 electric field measurements at high latitudes 
from a large nrimber of orbits. In this study the polar cap was 
taken to be the region of dawn- dusk electric fields, and the 
auroral zone the adjacent lower latitude region of poleward and 
equatorward electric fields in the evening and morning sectors 
respectively. While the measurements showed a large variability at 
the daytime boundary location, it was found that the overall high 
latitude electric field pattern is fairly consistent and does not 
fluctuate on a large scale during the incidence of substorms. 

A transition in the electric field from the evening to the morning 
pattern was foTind to occur in the premidnight 'sectctr -near .'2200 
hours geomagnetic local time. 

Meriwether et al. (1973) observed the motion of both ions and 
neutrals at high magnetic latitudes (> 65°) from electrically neutral 
luminous cloud and barium ion cloud pairs. These releases were 
made at various locations in the polar cap, the auroral belt and 
the transition region and the neutral cloud motions were determined 
using the photographic triangulation technique. The resultant height 
averaged wind vectors were plotted in geomagnetic latitude and 
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geomagnetic local time. The location of the auroral zone and the 
polar cap boundaries, as defined by Heppner (1972) was found to 
vary depending on convection conditions at the time of flight. As 
such, the choice of latitude and local time coordinates for certain 
of the wind vectors was somewhat a matter of discretion. Figure 1 
shows the relative positions of the neutral wind vectors and the 
auroral zone as determined statistically by Heppner (1972). 

Meriwether et al. observed a significant spatial rotation of 
the wind in the transition region between the auroral zone and the 
polar cap and m addition, a transition from westward to eastward 
convection near local midnight. With the exception of the early 
morning hours, it was found that the observed wind velocities 
could be modeled fairly well in terms of ion drag forces. 

A quantitative analysis of the variation of the Joule heating 
rate at various latitudes was performed by Ching and Chiu (1973) 
using measurements of ionospheric electric fields on a global basis 
(Haerendel, 1970; Mozer, 1971; Cauffman and Gurnell, 1971). This 
study showed the variation of Joule heating rates over a solar 
cycle and included both the altitude and latitude dependence as well. 
The calculated Joule heating rate at high latitudes was shown to be 
similar in magnitude and height distribution to the global solar EUV 
heating rate. Joule heating was thus thought to be a major 
contributor to the maintenance of the nighttime thermosphere. 

Using the incoherent scatter radar at Chatanika, Alaska 
(geomagnetic and geographic latitudes, 65.1°N, invariant latitude 
-A. = 65.5 N)^ Brekke and Doupnik (1973) have derived neutral 
winds for the E- region from measurements of ion drift motions for 
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both magaetically quiet and moderately disturbed days (Kp^ 3), The 
results represented the height range from 110 to 170 km and were 
modeled asstiming a day to night motion from the solar tide, as 
well as the ion drag force in the polar cap and auroral zone. The 
neutral atmosphere employed in the model was taken from Banks 
and Kockarts (1973). In summary, their observations indicated; 

(i) a strong poleward flow of the neutral winds near local 
noon and an equatorward motion in the morning (0700 
Alaskan Standard Time) and evening (2000 Alaskan 
Standard Time) hours. 

(ii) a clear correlation between electric field directions 
(deduced from ion motions near 170 km where ion- 
neutral coupling is minimal) and neutral wind directions, 

(iii) the importance of knowing the observers location relative 
to the auroral oval in making sense of observed wind 
vectors . 

Nagy et al. (1974) performed simultaneous neutral wind 
measurements from Fabry- Perot interferometer studies of the 
63 00 ^ atomic oxygen line at Ester Dome, Alaska and ionization 
drift measurements from the Chatanika radar station. Their results 
represented motions near 200 km with a spatial resolution of ion 
and neutral motions of about 150 and 400 km respectively in the 
horizontal -direction. In geomagnetic latitude and local time 
coordinates, the observations fell within the general region of 
convective motion from the nighttime to the daytime region, as 
determined statistically by Heppner (1972). The results conformed 



19 


to the general pattern observed in this region by Meriwether et al. 
(1973) and show a generally westward motion of ionization and the 
zonal component of the neutral wind in the pre-magnetic midnight 
sector. In the post midnight sector, the neutral gas did not, in 
general, follow the reversal of direction to the east as exhibited 
by the ion motion, implying either some additional unexplained 
force on the neutrals or a decrease in the ion drag force, 
presumably the effect of recombination. In order to explain 
occasional departures from the above general pattern, the authors 
suggested a combination of effects, including diurnal pressure 
gradients, energy input from auroral activity, J-oule heating and 
energy deposition from fast particles. 

Using a tmiform neutral density distribution and assmning 
height independent electron and ion densities, Maeda (1976[) ha^ 
solved the equations of motion for both ions and neutrals in the 
polar cap and auroral zone in the height range from 100 to 200 km. 
The electric field configuration used in the analysis was similar in 
magnitude and latitudinal variation to that observed by OGO-6 
(Gurnett, 1972). The model results were represented by polar plots 
of ion and neutral gas motion at altitudes of 110 and 200 km. 

Better agreement between the model and observations was 
found in the motion of the ions than that of the neutrals. The 
calculated neutral wind pattern, however, was found to depend 
critically on the degree of coupling between ions and neutrals and 
hence, on the assumptions made concerning the ion densities. In 
addition, certain of the observations used for comparative purposes 
were obtained on fairly disturbed days, whereas the density 
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distributions assumed in the model corresponded to a relatively 
quiet state of the polar ionosphere, 

1.2,4 Dynamical Heating at Thermospheric Heights 

In addition to solar energy input, a major source of heating 
at thermospheric heights has long been thought to be due^ to 
dynamical heating from atmospheric tidal motions of various scales. 

One of the early discussions pertaining to these tidal motions 
was given by Hines (I960), Hines asserted that a variety of 
irregular motions observed in the upper atmosphere can be 
explained in terms of upward propagating gravity waves generated 
low in the atmosphere and deriving their energy from either tidal 
oscillations and their harmonics or the low altitude wind system. 
Hines discusses gravity waves in terms of propagating adiabatic 
waves and cites direct evidence for the existence of these waves as 
shown by the results of early scattering experiments and height 
variations of the displacement of meteor trails. 

It was easy to account for these effects in terms of a wave 
theory since the wave interpretation accounted for observed large 
horizontal velocity components compared with vertical ones and, 
in addition, the large horizontal to vertical wavelength ratio of 
these disturbances. 

According to the theory developed by Hines (I960), the 
general pattern of these disturbances is revealed by a cyclic 
variation of wind direction with increasing height. Phase propaga- 
tion is in a downward direction and energy propagation generally 
upward, the energy of a given mode being removed by damping 



Z1 


effects as height increases. This effect is greater for modes with 
shorter periods and thus the longer period modes dump their 
energy at preferentially higher altitudes. 

Hines (1965) reviewed some aspects of his earlier theory in 
the light of observational data obtained after I960. From an 
analysis of mid latitude wind profiles obtained from rocket vapor 
trails, Kochanski (1964) found a relationship between the density p, 
the height z and the wave induced horizontal wind speed given by 

2 

pU = Const. • exp(-9.26z/z ) z = 140 km (1.2) 

X o o 

2 

This exponential decline in the product p tended to confirm the 
losses discussed by Hines (I960) and in addition, dominant modes 
observed both by Kochanski (1964) and by Zimmerman (1964). This 
decline also gave excellent agreement with the theoretical cutoff of 
smaller wavelength modes by viscous damping. 

Hines stated that gravity waves may profoundly affect the 
structure of the entire thermosphere and may in fact be the source 
of the rapid increase in temperature with height in the region 
from 100 to 110 km. As a source of heating at F-region heights, 
particularly at times of high geomagnetic activity, Hines proposed 
that the energy of gravity waves at high geomagnetic latitudes 
propagates horizontally over large distances to the low and mid 
latitude thermosphere, and is the source of heating observed by 
satellite measurements (Jacchia, 1964). In discussing the energy 
input to the E- region from the semi-diurnal tidal flux, Hines (1965) 
suggested that the semi-diurnal tide may be created in- situ and 
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that it may, in fact, be the diurnal tide and the prevailing winds 
that are responsible for electrodynamic heating at higher altitudes. 
This has been discussed more recently by Salah and Evans (1973). 

Direct confirmation of the presence of wavelike density 
variations in the thermosphere was provided by Newton et al. (1969) 
using density gauges aboard the Explorer 32 satellite. Data was 
obtained in the altitude range from 286 to 510 km with a spatial 
resolution along the orbital path of about 16 km. Their results 
indicated a general north- south phase propagation with maximum 
horizontal wavelengths of several hundred km, in agreement with 
earlier analyses by Thome (1964) and Georges (1968). These 
observations were best explained in terms of gravity waves 
propagating predominantly horizontally in a north to south direction 
with horizontal wavelengths between 130 km and 520 km. 

Observations of the F- region from incoherent scatter at 
St. Santin - Nanjay during magnetically disturbed conditions have 
been performed by Testud (1970). These results revealed an 
enhancement of oscillations in the wind and temperature structure 
following magnetic events, with a time delay of about 1 hour. These 
disturbances have been interpreted in terms of gravity waves, and 
simultaneous observations by the European network of ionospheric 
observatories showed that the disturbances propagate southward 
(Testud and Vasseur, 1969). They are thus thought to be of auroral 
origin, in agreement with the theory of Hines (1965) and Newton et 
al. (1969). A theoretical analysis by Testud once again supported 
the view that large scale travelling ionospheric disturbances (TID) 
are generated at high latitudes and transport energy to mid and 
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low latitudes . 

Klostermeyer (1973) performed a theoretical study yielding an 

estimate of the height distribution of reversible heat input from 

gravity waves and the resulting temperature increase of the neutral 

gas. The calculations were performed using parameters typical of 

both quiet and disturbed conditions. In the quiet condition set, the 

author employed the 1966 U.S. Standard Atmosphere Supplements 

assuming an exospheric temperature of 900 K. The ion number 

density distribution was represented by an n- Chapman layer with a 

11 -3 

maximum density of 5x10 m at 300 km. The height distribution 
of heat input from a single wave mode was found to be similar to 
that of solar EUV heating, as presented in a theoretical study by 
Volland (1969) and resulted in an increase of about 40 K in the 
exospheric temperature over a period of about 2 hours. 

The calculations for increased geomagnetic activity employed 
a perturbed density distribution and an expression relating the 
exospheric temperature to the Kp index (Jacchia et al. , 1967). The 
results showed good agreement with observed temperatures for Kp 
levels ^ 5, and best results occured when a time delay between 
disturbance and effect of about 6 hours was assumed, implying 
speeds of horizontal propagation of the waves comparable to the 
local sound speed. 

A study of the effects of viscous and thermal dissipation on 
horizontally ducted acoustic-gravity waves was given by Francis 
(1973a). An important difference between the analysis by Francis 
and earlier ones was the use of a more realistic sound speed 
profile in the thermosphere. The author interpreted traveling 
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ionospheric disturbances (TID) as upper atmosphere gravity modes 
detected by virtue of their effect on the F- region electron density. 
Included as a subset of this broader classification are those waves 
generated in the auroral zone during magnetic substorms. The 
results confirmed that large scale TID with speeds greater than 
300 m sec”^ are ducted horizontally for large distances with 
varying degrees of attenuation, depending on the vertical component 
of the energy flux. In addition, the author discussed the tempera- 
ture dependence of the horizontal phase velocity in the upper 
thermosphere. The velocity was found to vary by roughly a factor 
of 2 between the extremes of temperature encountered and was 
larger for a given mode during daylight than at night. 

In a later paper, Francis (1973b) suggested that medium scale 
TID are manifestations of freely propagating gravity waves rather 
than ducted modes, used to explain large scale TID, and that 
the direction of propagation is constrained to follow the curvature 
of the earth as a result of refraction by the earth’s gravitational 
field. 

In addition to gravity waves with periods above the Brunt 
Vaisala cutoff and below the viscous dissipation limit, there exist 
atmospheric oscillations related to the solar tide. These periodic 
oscillations are also revealed in temperature and density profiles 
as well as in the winds, and involve a fundamental period of one 
solar day and the various harmonics. In tidal theory these 
oscillations are represented by an expansion in spherical harmonics 
over the indices 1, m and n. It is customary to designate a given 
term by where the indices 1 and m classify the term as to 
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wave family and the n index further subdivides it into wave type. 
The index 1 specifies the period of the wave, given by 1 ^ of a 
mean solar day, and 1 and n together (1, n) are normally used to 
specify the tidal mode under consideration. 

Using incoherent scatter data obtained at Millstone Hill 
(42.6°N) and St. Santin, France (44.6°N), Salah and Wand (1974) 
have studied the average behavior and variability of tidal 
oscillations in the lower mid latitude thermosphere. Good agree- 
ment was obtained between the results from the two stations. 
Analysis of the data revealed that the dominant tidal component at 
mid latitudes in the lower thermosphere is the semi-diurnal (2,4) 
mode. 

The semi-diurnal oscillations can be resolved into two semi- 
diurnal tidal modes, the (2,2) and the (2,4) modes. The results of 
Salah and Wand indicated an attenuation of the tidal (2,4) mode in 
contrast with theory (Butler and Small, 1963; Richmond, 1971; 
Lindzen and Hong, 1974), yet the contribution of the (2,2) mode at 
heights between 100 and 125 km is significantly less than predicted, 
allowing the (2,4) mode to dominate in this region. Additional 
disagreement between theory and observations was found in the 
height of maximum amplitude of the two modes and in the damping 
rates above these heights. 
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1. 3 Specific Sfcatemenl: of the Problem 

It IS the intent of this study to; 

1) provide estimates of the energy transported into or out of 
the nighttime thermosphere via conduction, radiation and 
the low and mid latitude thermospheric wind system. 

2) estimate the global energy input from the polar cap to the 
mid latitude nighttime thermosphere via the high latitude 
wind system. 

3) ascertain the need for additional thermospheric energy 
sources based on a comparison of the energy input and 
loss terms. 



CHAPTER II. 


METHOD OF ANALYSIS 

2.1 Analytic Form of the Contributions to the Energy Balance 

For a composite of' several neutral constituents i, the thermal 
balance of the system under steady state conditions requires the 
solution of the energy equation. The form used here is that given 
by Banks and Kockarts (1973, equation 14.12). 
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where the heat flow vector E is 
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In addition, the continuity equation for the total density is given by 
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In the above expressions 

v^ is the mass average or bulk velocity 

n^^ , m. , and are, respectively, the concentration, mass 

I 

and diffusion velocity of the i th constituent 
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n is the total concentration of neutral constituents 
T is the temperature 
X. is the thermal conductivity 
p is the total hydrostatic pressure 
k is the Boltzmann constant 

F. are external forces acting on the i'' th , constituent 
N is the average number of degrees of freedom for the 
constituents, defined by 


N 


2N.n. 
■1 ^ 1 
Sn. 
i 1 


n 



1 


(2.4) 


The work term 



in equation (2.1) 


can be 


rewritten as 


pV* v^ = V. (pT^) - T^-Vp 


(2.5) 


Ass\iming hydrostatic equilibrium, equation (2.5) becomes 


pV*v^ = V.(nkTv^) - v^. (n m g ) 

Writing 


( 2 . 6 ) 


g = - (2.7) 

where <p is the gravitational potential, equation (2.6) becomes 

p^^.Vo = V.(nkT'v^) + 'v^.(nm'^^P) (2.8) 
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and expanding gives 


= V*(nkTv^) 4- • (n m v^ <jP ) -■ ^PV.(nmv^) (2.9) 

From equation (2.3) and for steady state conditions, 

"V.(nmT^) = 0 (2.10) 


Thus, equation (2.6) can be rewritten as 



V.(nkTv^) + V.(nmv^<P) 


( 2 . 11 ) 


In performing the heat balance calculations, the zero of 
potential energy has been taken to lie on the lower system 
boundary at 120 km, as shown in Figure 2. For a flat earth 
approximation, the gravitational potential energy term <p in 
equations (2.8) through (2,11) is given by g(h-h^). This approxi- 
mation is within 6 percent of the true value at h= 300 km and 
within 2 percent at h= 150 km. Making this substitution and 
integrating equation (2.1) over the entire system volume, one gets 
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where Q 

R 


IS the total radiative energy loss rate 


m is the mean molecular mass 

W is the difference between the total observed losses and 
the energy sources shown 


The term 


I 


n. F. 



represents the work done on the type i 


i 

particles by the vector sum F^^ of the external forces. The only 
external forces that contribute to this term are frictional forces 
arising from the relative motion between the neutral gas and ionic 
species present in the gas. This interaction is commonly labeled 
ion drag. 

It is among the purposes of 'this study to establish the need 
for additional energy source or loss terms based on estimates of 
the conduction, radiation and the thermospheric wind system under 
steady state conditions. These additional terms can be accounted 
for in terms of either ion- drag ‘or the quantity W on the right hand 
side of equation (2.12). The thermal balance, then, requires an 
estimate of each of the following established energy input and loss 
terms, expressed as fluxes across the system boundaries. 


Energy input terms: 


(i) thermal kinetic energy of the major constituents 

- ^ nkT^^ + ^(^ + l)n.kT V. ] • d’S (2.13) 

‘'s' i 
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(li) gravitational potential energy with respect to the lower 
boundary at 120 km 


^ n m g (h - h^) ^ 


(2.14) 


Energy loss terms: 


(iii) downward conduction {Qq) a-t 120 km due to the 

vertical temperature gradient 8T/9z at that height 
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(iv) infrared radiative cooling 


Q 
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The energy input terms arise mainly from the transport of 
neutral gas into the nighttime region via the horizontal thermo- 
spheric wind system. The problem of the energy flow, then, 
requires a simultaneous mass balance, and this entails the 
solution of the mass continuity equation (2.3). For a steady state 
process, the integration of equation (2.3) over the volume of the 
system yields 
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Thus, a mass flux into or out of the system from the thermo- 
spheric winds gives rise to a compensating flux at the 120 km 
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boundary. This boundary flow, in turn, results in an accompanying 
energy flux across the lower boundary which must also be " 
considered in order to arrive at a net energy gain to the system. 

The transport of energy into the thermosphere via gravity 
wave propagation and in addition particle precipitation and Joule 
heating processes at high latitudes has been well established 
(Hines, 1965,1969; Roper, 1966; Justus, 1967). These processes 
constitute independent energy sources in addition to those considered 
above. Recent estimates (Klostermeyer , 1973) show that thermo- 
spheric heating from these sources is comparable in size with the 
above mentioned processes. The difference between the total energy 
input and loss rates given by items (i) through (iv), then, should 
provide an independent estimate of the thermospheric heating from 
the less well understood processes in the nighttime region. 

2. 2 The Coordinate System Employed 

The choice of suitable boundaries for the nighttime thermo- 
sphere is influenced by a number of factors. Estimates of the 
fluxes across the boundaries require a distribution of experimental 
data over at least part of the surface, and hence the choice must 
correspond to the locations of the largest body of available 
measurements . 

The geometry involved in performing the energy balance of 
the nighttime thermosphere for equinoctial conditions is shown in 
simplified form xn Figure 2. The system is defined as the anti- 
solar region boxinded below by a spherical surface 120 km from the 
surface of the earth, and in front by the morning and evening 
terminators. At equinox, the polar axis and the terminator lie in a 
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Figure 2 


Geometry employed in performing the energy balance of the nighttime thermosphere 
at equinox. 
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common plane and the corotating component of the total horizontal 
wind is normal to the terminator at each point. At other than 
equinoctial conditions, the polar axis is inclined at an angle to the 
plane containing the terminators, this angle being a maximum of 
23.5 degrees at s\immer and winter solstice. The upper boundary 


is typically taken as 


the height z where the quantity 


UJ 


kTdz' 


has 


z 

diminished to a small fraction (<.01) of the value at 120 km. This 
insures that the contributions to the energy balance from heights 
above this level are negligible in comparison. In practically all 
cases this requirement has been met at an altitude of 300 km. 

As shown in Figure 2, the nighttime thermosphere considered here 
is taken to be fixed with respect to a coordinate system having the 
earth- siin line as one axis. The motion of the neutral gas across 
the terminator boundaries as viewed in this fixed frame is the 
vector sum 


^fixed 


V . + u X r 

moving 


(2.17) 


where v . is the velocity of the neutral wind observed in the 
moving •' 

rotating frame. These relative velocities are 
measured or deduced using rocket vapor trail and 
incoherent scatter techniques, 
w is the angular velocity of the earth's rotation 




fixed 


is the velocity of the neutral gas in the fixed frame 


of reference. This is the velocity used in the 
computation of the horizontal energy and mass fluxes 
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and is related to v and V. discussed earlier by 

o 1 


fixed . 


V + V. 
o 1 


(2.18) 


for the i'th constituent. 

r IS the radius vector from the center of the earth to the 
location of the observed air parcel. 


The polar regions have been excluded from the nighttime 
system defined here since the energy considerations are completely 
different as discussed by Nisbet and Glenar (1977). While the high 
latitude boundaries are defined by the cutaway portions at the top 
and bottom of Figure 2, the polar winds are ordered in geomag- 
netic local time and latitude. It has been assumed for the present 
analysis that the average energy transported by the winds across 
the 70 degree geographic latitude into the nighttime thermosphere 
IS equal to the energy transported by the winds across the 70 
degree geomagnetic latitude over the corresponding range of 
geomagnetic local times. These south directed energy and mass 
fluxes constitute the energy input to the system from auroral 
processes and are treated in detail in Section 3.2. 

An additional complication arises at other than equinoctial 
conditions when a larger or smaller fraction of the polar cap 
area is contained within the region of interest. In such a case, 
the geometry shown in Figure 2 no longer applies. As discussed 
in a later section of this study, the small number of wind 
measurements made at high geomagnetic latitudes are not sufficient 
to show seasonal changes in the high latitude wind pattern. For 
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this reason, and since the terminator is always within a few 
degrees of the polar cap, the high latitude energy source is in 
each case taken to be the energy outflow from the cap between 18 
and 6 hours geomagnetic local times. 

Figure 3 shows pictorially the contributions to the energy 
balance described in Section 2,1. The treatment of the energy input 
from thermospheric winds has been divided into contributions from 
the low and mid latitude wind system and contributions from the 
polar regions. For this reason, the manner of presentation of 
contributing terms differs slightly from that listed in Section 2.1. 
Each quantity in the figure related to the horizontal wind contains 
both the thermal and gravitational potential energy contributions 
to the total flux. The contributing terms in Figure 3 may be 
defined as follows; 

Qp = net energy gain to the nighttime system from the 
polar thermospheric wind system 

(^p = total mass flow entering the system through the 

polar boundary, induced' by the polar winds 

Qp = net energy gain to the system from the mid and low 
latitude thermospheric winds 

C|>P = total mass flow entering the system, induced by the 
mid and low latitude winds 

Some additional mention should be made regarding the 
definition of net energy gain as it is used here. As shown in 
Figure 3, the net energy gain represents only the portion of the 
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Figure 3 Contributions (symbolic) to the energy balance from the energy input and loss 
processes. 
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energy conveyed by the horizontal winds across the system 
boundaries which may be considered as energy input in performing 
the heat balance. By the mass continuity equation (2.16), the mass 
inflow through the terminator or polar boundary is assumed to 
leave the system at the 120 km boundary with some suitable 
estimate for the 120 Ion temperature. The energy leaving in this 
manner must be subtracted to arrive at a net energy gain. Each of 
the terms above is discussed in detail in a later section of this 
study. 



CHAPTER ni 


CONTRIBUTIONS TO THE ENERGY BALANCE FROM 
THE THERMOSPHERIC WINDS 


3.1 Energy Transport Via the Low and Mid Latitude Wind System 
As shown in Figure Z, the low and mid latitude thermospheric 
winds convey energy into or out of the nighttime system at the 
terminator boundaries and it is the component of the horizontal 


wind normal to the morning and evening terminators that 

contributes to this exchange. The portion of the energy surface 

integral over the morning and evening terminators » expressed 

in polar coordinates, is given by equation (3.1). The diffusion 

velocities V. have been found to be small compared to the total 

2 2 1/2 

horizontal wind velocity (v + V. ) and the contribution from V. 

•' O I 1 

in equations (2.12) and (2.13) has thus been omitted here, 
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(3.1) 


where is the component of the total horizontal wind velocity v^ 

normal to the terminator (positive into the system) 

R^ is the radius of the earth (6370 km) 

S,,,-. indicates the sum over both terminators 
ME 


The second term on the right hand side of equation (3.1) arises 
from a consideration of the mass continuity equation (2.16) and 
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represents the thermal energy lost at the lower boundary. Again, 
is the mass flux delivered to the system by the horizontal winds. 

R^+300 70° 

®MeI r C (nim.vj dvrdr (3.E) 

i R +120 -70° 
o 

It is clear from equations (3.1) and (3.2) that an estimate of 
the potential and thermal kinetic energy flow across the terminator 
boundaries requires a knowledge of the height and latitudinal 
dependence of the neutral temperature, the major constituent number 
densities and the thermospheric wind velocities. 

3.1.1 Sources of Wind Data; Vapor Trail Measurements 

Rocket vapor trail measurements provide information on the 
magnitude and direction of the horizontal winds at heights up to 
nearly 200 km. The results of these flights are particularly 
valuable for the purposes of this study since the technique involves 
a direct measurement of the wind velocity and since a large 
number of the rocket flights are made at or near the terminator 
location so that trails illuminated by the sun can be observed 
against a dark background. 

In this analysis, data has been examined from a number of 
morning and evening pairs of sodium and lithium vapor trail 
measurements made at mid latitudes between the years 1961 and 
1972 . A summary of these flights and associated observations is 
given in Table 1. In grouping the results according to season, 
data obtained within thirty days of either equinox was considered 
representative of that equinox, while the remaining firings were 



Table la. Mxd Latitude Vapor Trail Measurements 
Near Fall and Spring Equinox 


Date 

Mo/Day /Yr 

Day 

Time 
of Day 

Launch Site 

t 

K 

P 

F 

10.7 

Altitude Range 
(km) 

Ref. 

9/7/61 

250 

PM 

Sardinia, Italy 

0.6 

115 

90-200 

A 

9/8/61 

251 

AM 

39.6°N, 9.5°E 

0.8 

117 

90-140 


9/16/61 

259 

PM 

Wallops Island, Va. 

2.0 

133 

100-150 

A 

9/17/61 

260 

AM 

37.8°N, 75.5°W 

3.4 

124 

100-170 


10/7/64 

280 

PM 

Wallops Island, Va. 

1,6 

74 

90-150 

A 

10/8/64 

281 

AM 


2.8 

77 

85-160 


10/16/69 

289 

AM 

Woomera, Australia 

0.8 

116 

95-200 

B 

10/17/69 

290 

PM 


1.5 

123 

95-240 


4/19/61 

109 

PM 

Sardinia, Italy 

1.4 

105 

85-185 

A 

4/20/61 

110 

AM 

39.6°N, 9.5°E 

2.5 

.03 

110-190 


4/19/61 

109 

AM 

Wallops Island, Va, 

2.0 

105 

90-160 

A 

4/20/61 

110 

PM 


1.0 

103 

80-160 


4/21/61 

111 

AM 


0.2 

104 

85-160 


3/1/62 

60 

PM 

Wallops Island, Va. 

2.1 

121 

70-125 

A 

3/2/62 

61 

AM 


1.9 

112 

65-125 


4/12/72 

102 

PM 

Elgin AFB, Fla. 

1.4 

128 

104-168 

C 

4/13/72 

103 

AM 

30.4°N, 86.7°W 

2.9 

128 

123-174 


4/14/72 

104 

PM 


2,5 

124 

82-132 


4/14/72 

104 

PM 


2.5 

124 

183-194 



References: A. Bedinger (1966), B. Lloyd et al. (1972), C. Roper and Edwards (1974), D. Justus and 

Edwards (1965), E. Rees et al. (1972), F. Montgomery et al. (1968). 

^ Average of three values preceding measurement. 



Table lb. Mid Latitude Vapor Trail Measurements 
Near Summer Solstice 


Date 

Mo/Day/Yr 

Day 

Time 
of Day 

Launch Site 

+ 

K 

P 

F 

U0.7 

Altitude Range 
(km) 

Ref. 

6/5/62 

156 

AM 

Hammaguir, Algeria 

2.8 

85 

100-150 

A 

6/5/62 

156 

PM 

31°N, 3°W 

1.6 

85 

105-140 


6/5/62 

156 

AM 

Reggan, Algeria 

2.8 

85 

95-140 

A 

6/5/62 

156 

PM 

26.70N, 0° 

1.6 

85 

95-120 


5/20/63 

140 

PM 

Sardinia, Italy 

0,5 

91 

95-195 

A 

5/21/63 

141 

AM 

39.6°N, 9.5°E 

0,5 

88 

95-160 


5/21/63 

141 

PM 


1.6 

88 

85-180 


5/17/63 

137 

PM 

Elgin AFB, Fla. 

0.6 

100 

95-170 

A 

5/18/63 

138 

AM 

30.4°N, 86.7°W 

0,4 

98 

100-160 


5/17/63 

137 

PM 

Elgin AFB, Fla. 

0.6 

100 

71-170 

D 

5/18/63 

138 

AM 


0,4 

98 

99-105 


5/21/63 

141 

PM 




69-140 


5/31/68 

151 

AM 

Woomera, Australia 


145 

95-240 

E 

5/31/68 

151 

PM 

31°S, 137°E 


145 

95-240 



References: A. Bedinger (1966), B. Lloyd et al. (1972), C. Roper and Edwards (1974), D. Justus and 

Edwards (1965), E, Rees et al. (1972), F. Montgomery et al. (1968). 




Average of three values preceding measurement. 



Table Ic. Mid Latitude Vapor Trail Measurements 
Near Winter Solstice 


Date 

Mo/Day/Yr 

Day 

Time 
of Day 

Launch Site 

K ' 
P 

F 

10.7 

Altitude Range 
(km) 

Ref. 

11/27/62 

331 

PM 

Chamical, Argentina 

1.7 

75 

105-170 

A 

11/30/62 

334 

AM 

30.5°S, 66^W 

3.5 

77 

100-170 


11/30/64 

334 

PM 

Sonmiani, Pakistan 

3.4 

74 

120-155 

A 

12/1/64 

335 

AM 

26°N, 67°E 

1.9 

76 

100-205 


11/11/64 

315 

PM 

Ship at 34.4°N, 74,7°W 

0.8 

72 

80-155 

A 

11/12/64 

316 

AM 

34.4°N, 76.3°W 

2,4 

72 

90-155 


11/17/65 

321 

PM 

Elgin APB, Fla, 

0.5 

74 

97-155 

F 

11/18/65 

322 

AM 

30.4°N, 86.7°W 

1,5 

75 

94-152 


11/18/65 

322 

AM 


1.5 

75 

89-171 


11/18/65 

322 

AM 


1.5 

75 

78-138 


1/16/67 

16 

AM 

Elgin AFB, Fla (Upleg) 

2.7 

124 

86-158 

F 

1/16/67 

16 

AM 

(Downleg) 

2.7 

124 

100-156 


1/16/67 

16 

PM 

(Upleg) 

1,7 

124 

88-162 


1/16/67 

16 

PM 

(Downleg) 

1.7 

124 

92-102 



References: A, Bedinger (1966), B. Lloyd et al. (1972), C. Roper and Edwards (1974), D. Justus and 

Edwards (1965), E. Rees et al, (1972), F, Montgomery et al» C1968). 


t 


Average of three values preceding measurement. 
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considered representative of either summer or winter solstice. 

With the exception of the two pairs of Australian rocket firings 
shown in Tables lb and Ic, all of the morning - evening rocket 
flights were performed in northern latitudes. For computation 
purposes the average wind profiles in Figure 4 are thus considered 
to represent only the northern hemisphere. 

For each date involved, the wind flow into the system is 
evaluated by first considering the angle between the east-west 
direction and the terminator and then evaluating the component of 
the wind normal to the terminator. By averaging the results 
obtained at each terminator in increments of 5 km, irregularities 
in the wind system from propagating waves of various scales has 
been largely removed and a consistent trend in the variation of wind 
propagation with height is revealed (Figure 4), 

The effects of changing magnetic activity have been taken into 
account by monitoring the magnitude and daily range of the Kp 
index for the date of each firing. The observed Kp values are all 
characteristic of low or medium levels of magnetic activity and in 
no case is there a wide daily range in Kp indicating a sudden 
increase in magnetic activity between the morning and evening 
firings. Thus, the contribution from this variable to the daily wind 
variations is thought to be small. 

Another of the uncertainties involved in these calculations is 
the effect of variations in solar activity as indicated by the F,„ „ 

J. U • I 

index, on the energy transported into the nighttime region and 
hence on the observed zonal fluxes at the terminators. A density 



ALTITUDE (KM) 


200 — 


( 8 ) ( 6 ) 



a) MORNING PROFILE (v^) 
(m sec“’} 


(lOKB) 



Figure 4 Height variation of the average eastward neutral wind at twilight, from mid- latitude 
vapor trail measurements . The numbers in parentheses are the number of flights used 
xn obtaining the respective profile. 
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weighted average velocity above 120 km was defined by 



(3.3a) 


:(3*.3b) 


where p , H and T are the density, mean scale height and 
temperature at height z. Here, Vj_ is defined such that a positive 
value of y^ indicates wind flow from the dayside to the nights ide 
hemisphere. This quantity should provide an indication of the 
direction of height integrated mass and energy flow above 120 km. 

To check the effect of solar activity on the values of V|_, a 
correlation study was made between the values of y^ for each flight 
and the ^ value at the time of each flight. Representative 

temperature distributions were taken from the COSPAR International 
Reference Atmosphere (CIRA, 1972). No significant correlation 
between the values of Vj^ and the Fj^q ^ index was observed, and 
since the range in observed values of Fj^q ^ was not excessive 
(72- 145), no restrictions were placed on values of the solar index 
in selecting applicable data. 

The profiles obtained by averaging wind data for the same 
season are shown in Figure 4. Figure 4c has been included to 
illustrate the height variation of mass and energy transport from 
the dayside to the nighttime region. At heights above 140 km, these 
profiles invariably show, as expected, an anti - solar flow, 
presumably the result of dayside solar heating. With the exception 
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of the summer solstice results, the direction of the net wind 
curiously reverses in the region near 135 km, and the density 
weighted average velocity vj^ is found to be -22 and -20 m sec ^ 
for equinox and winter solstice periods. For comparative purposes, 
a similar analysis was performed on characteristic summer and 
winter wind profiles deduced by Kochanski (1964), and extending to 
140 km. These were based on 25 mid latitude morning and evening 
twilight rocket vapor trail measurements. 

A negative value for was obtained using Kochanski' s 
summer data, but will not be discussed since the data is not 
completely independent of that .used in this study. All of the 
morning and evening flights grouped into Kochanski' s winter 
classification are different from those used here, however, and 
application of equation (3.3b) to these results shows a value for 
of -15 m sec \ in reasonable agreement with the winter solstice 
value stated above. 

At heights above 140 km. Figure 4 shows a westward flowing 
wind at the morning terminator and an eastward flowing wind at the 
evening terminator with an accompanying day to night energy flow, 
in agreement with tidal theory. A significant change in wind 
direction is evident at lower altitudes, however, particularly in the 
equinoctial results. The effect of this is to produce a complete 
turn-around in the direction of the height integrated mass and 
energy flow at the terminator as the lower integration limit is 
moved from 135 down to 120 km. 
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3.1.2 Analysis oi Temperature and Density Data from 
Arecibo and St. Sanfcin 

In order to calculate the energy transport across the morning 
and evening terminators, it is necessary to make use of the wind 
results from the previous section. In addition, estimates of the 
densities and temperatures at both terminators and at the altitudes 
of interest are required. An additional source of data, then, are 
incoherent scatter measurements of densities and temperatures at 
the terminators as provided by the Arecibo {18.3°N) and the 
St. Santin (44.6°N) incoherent scatter facilities. 

The Arecibo facility provides measurements of the total 
collision frequency at selected heights in the altitude range from 
90 to 110 km which have been scaled to provide the total collision 
frequency and hence the total density at 120 km. In order to 
estimate the atomic oxygen number density at 120 km, it was 
assumed that the CIRA, 1972 model represents the atomic oxygen 
density correctly at 320 km and that diffusive equilibrium allows 
the calculation of the 120 km density with the aid of the tempera- 
ture profile discussed below. The molecular number density was 
divided between and N^, in the ratio 1:7, again following CIRA, 

The temperature profile at the altitudes of importance is 
deduced using observed ion temperatures at heights between 118 and 
260 km where it is assumed that the neutral temperatures do not 
depart significantly from that of the ions. An additional neutral 
temperature at 320 km is obtained according to the procedure of 
Swartz and Nisbet (1971) from the CIRA, 1972 densities and the 
observed ion and electron temperatures. Bates (1959) represents the 
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height distribution of the neutral temperature by 

= 'Tco “ exp'(.s«^>) (3.4) 

where is the neutral temperature at height z 

IS the exospheric temperature 

is the neutral temperature at some reference altitude 

z j here taken to be 120 km 
o 

s is the shape parameter 
<p is the geopotential altitude 
(z - z ) (R. + z ) 

<p = z - z for z « R (3.5) 

R + z ° 

Equation (3.4) has been solved using the least squares technique 
to obtain the best fit values of the temperature parameters 

and s. The uncertainty in is not of any consequence here, 
since only temperatures well below 200 km are important in the 
computation of energy fluxes height integrated above 120 km. 

Table A.l summarizes the temperatures and densities derived 
in the above manner using data from 18 days observations at 
Arecibo between August 12, 1974 and November 17, 1976. The data 
used in this analysis represent low or medium levels of magnetic 
activity (Kp ^ 4.0). The dates were selected from a body of 41 days 
observations and represent those days where runs were made 
within 1.5 hours of both terminators and where the error in the 
observed collision frequencies was less than 50 percent of the 
stated value. By examining results of the first few observations in 
the morning and the last few in the evening of each day, one can 
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infer the temperatures and densities at the location of the 
terminators. Since each of the observation days listed has been 
classified according to season, together with the appropriate wind 
profile (Figure 4), the values should provide the information 
necessary to determine the low latitude contribution to the energy 
surface integral, expressed by equation (3.1), 

Figure 5 shows the calculated energy flux per unit terminator 
length at the time of each observation and the extrapolated flux 
across the terminators for each of the observation days considered. 
The term "net flux" is used to represent these values since a 
correction has been applied in each case for the lower boundary 
flux, as discussed at the beginning of this section. The base 
temperature used in applying this correction is the average of the 
120 km temperatures at the morning and evening |terminators for 
the appropriate observing day. The terminator fluxes have been 
determined using a linear least squares fit to the data. The error 
bars represent the \mcertainty in each determination which arises 
from both the error in the derived collision frequencies and the 
intrinsic variability with local time of the densities in this altitude 
range. Table 2 lists in siimmary form the terminator fluxes 
obtained for each of the 18 'observing days. The quantities AW and 
AR in the table represent the contributions to the daily net flux 
from the measured winds and the corotating winds respectively. 

As can be seen in Figure 5 and Table 2, both the magnitude 
and sign of the net energy flux vary between measurements. When 
grouped according to season, the average net fluxes were found to 
be 0.28, -0.62, -0.93 and -2.20X10^ W m”^ for summer solstice. 
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Table 2. Summary of Arecibo Results 


Net Flux Into 


Nighttime 


Date 

L.S.T. 

Season 

Thermosphere 
(N m-1) 

Daily Net Flux 
AW+AR (N m“^) 

Error 
(W m“^) 

Aug, 12, 1974 

5.66 

18.34 

Summer 

-l.lxioj 

4.8x10 

3.7x10^ 

5.8xl0'^ 

Aug. 14, 1974 

5.67 

18.33 

Summer 

-4.9x10^ 

5.1x10^ 

0.2x10'^ 

2.3x10^ 

Sept. 17, 1974 

5.95 

18.05 

Fall 

-5.0x10^ 

2.6x10^ 

-2 . 4x10^ 

1 . 6x10^ 

Sept. 18, 1974 

5.96 

18.04 

Fall 

-4 . 3xlot 
4.4x10^ 

0.1x10^ 

2.2x10^ 

April 15, 1975 

5.78 

18.22 

Spring 

-6.9xl0t 
4 . 4x10 

-2.5x10^ 

2.0x10^ 

April 16, 1975 

5.77 

18.23 

Spring 

-5. 4x10 t 
3.2x10 

-2.2x10^ 

3.2x10^ 

May 22, 1975 

5.53 

18.47 

Summer 

-3.4xlot 

4.5x10 

1.1x10^ 

1.4x10"^ 

June 10, 1975 

5.46 

18.54 

Summer 

-4.6x10^ 

3.2x10 

-1.4xl0'^ 

2.0x10^ 

June 11, 1975 

5.46 

18.54 

Summer 

4 

-3.8x10, 

2.9x10^ 

-0.9x10^ 

1.5x10 

Oct. 14, 1975 

6.18 

17.82 

Fall 

-5 .4xlot 
2 . 3x10 

-3.1x10^ 

2 . 8x10“^ 


AR 

0.22 

0.12 

0.19 

0.92 

0.60 

0.29 

1.84 

0.23 

0.50 


-4 


0.11 



Table 2. (Continued) 


Date 

L.S.T. 

Season 

Net Flux Into 
Nighttime 
Thermosphere 
(W m~^) 

Daily Net Flux 
AW+AR (W m"3.) 

Error 
(W m“^) 

1 ^^1 
1 ar| 

Jan. 20, 1976 

6.47 

Winter 

-6 . Ixiot 

/ 




17.53 


5 . 3x10 

-0.8x10^ 

1.3x10 

0.23 

Jan. 21, 1976 

6.46 

Winter 

-5.0x10^ 

A 

L 



17.54 


2.9x10 

-2,1x10 

1.4x10^ 

0.13 

Mar. 17, 1976 

6,02 

Spring 

-7,lxlof 

A 

4 



17.98 


5,2x10 

-1,9x10 

2.6x10 

0.29 

May 13, 1976 

5.57 

Sunimer 

-4.4x10^ 


A 



18.43 


3.4x10^ 

-1 . 0x10 

1.7x10^ 

0.20 

Sept. 22, 1976 

6.00 

Fall 

-2, 2x10 J 

/, 

4 



18.00 


4,4x10 

2.2x10 

1.3x10 

0.19 

Oct. 19, 1976 

6.23 

Fall 

-2.8x10^ 

A 

4 



17,77 


2.9x10 

0.1x10 

1 . 4x10 

6.00 

Nov. 16, 1976 

6.44 

Winter 

-3.2x10^ 

4 

4 



17.56 


3,1x10 

-0.1x10 

0.9x10 

0.57 

Nov. 17, 1976 

6.44 

Winter 

-2. 9x10 J 

A 

4 



17.56 


2.2x10 

-0.7x10 

1.2x10 

0.16 
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fall equinox, winter solstice and spring equinox respectively. These 

values suggest that the net day to night flux varies in a cyclic 

fashion over the course of a year and that at the latitude of Arecibo, 

energy tends to leave the nighttime region, particularly in the 

winter hemisphere. Because of the large uncertainties involved, 

however, the error in each of the above values is of the same 

order, or larger than, the value itself and any order which appears 

in the above seasonal fluxes must be viewed with caution. Using 

the standard deviation of the individual determinations for each 

season and the average values quoted above, the day to night energy 

flux at the latitude of Arecibo has been found to lie between - 

-1.6X10"^ and 2.2X10^ W m ^ for summer solstice, -2,.8X10'^ and 

1.5X10^ Wm ^ for fall equinox, -1. 8Xlo"^ and -0.1X10^ for winter 

4 4-1 

solstice and -2.5X10 and -1.9X10 W m for spring equinox. 

In order to examine the latitudinal dependence of the termin- 
ator energy fluxes, the wind results of Section 3.1.1 have also been 
used in connection with incoherent scatter results at mid latitudes. 
Data on the height distribution of electron and ion temperatures and 
the electron density at altitudes above 95 km, and the total 
collision frequency between 95 and 120 km are available from the 
St. Santin incoherent scatter facility (44.6°N). 

In the case of the St, Santin facility, each measurement does 
not furnish an entire height profile, but only data at a single 
height with a repetition period of about 1.2 hours for the E-region. 
Further, on a given observing day, E- region measurements are 
made from about sunrise to close to sunset, followed by F- region 
measurements made during the next 26 to 30 hours. To obtain data 
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at the terminator, the collision frequencies and temperatures for 
each height were plotted as a function of time and extrapolated 
visually to the time of sunrise and sunset. The total density at 120 
km was calculated from these values and used, together with the 
120 km temperature for the terminator flux calculations. Although 
the temperature profiles at the terminators were constructed from 
the F- region data from the following day, in no case was there a 
significant difference between the temperature parameters esti- 
mated for two successive mornings or evenings. 

Using this scheme, the St. Santin results were analyzed to 

obtain the temperature and density distribution at the morning and 

evening terminators during nine days in. 1969. and 1970.. These 

selected days are characterized by low and medium levels of 

magnetic activity and minimum variation in Kp values between the 

morning and evening hours. Table A.. '2- shows, the boundary conditions 

obtained from the St, Santin analysis and Table 3 lists the results 

obtained using the measured and corotating winds discussed 

previously and the listed boundary conditions. In contrast to the 

Arecibo findings, these results consistently indicate a net flux into 

the nighttime region. After grouping the results according to season, 

the energy flux into the nighttime thermosphere at St. Santin has 

been found to lie between 1.0X10 and 3.6X10 W m for summer 

solstice, 1.0X10^ and 3.8X10*^ Wm ^ for fall equinox and 1.0X10^ 

4 -1 

and 3.1X10 W m for winter solstice. 

It should be mentioned that in no case were measurements 
made at St. Santin before local sunrise or after local sunset. As 
a result, the derived temperatures and densities at the terminators 



Table 3. Summary of St 


Santin Results 


Net Flux Into 
Nighttime 





Thermosphere 

Daily Net Flux 

AW 

Date 

L.S.T. 

Season 

(W m“^) 

AW+AR (W m"l) 

AR 

Aug. 26, 1969 

5.30 

18.70 

Fall & 
Summer 

-2 . 7x10^ 
6.5x10 

A 

3 . 8x10 

0.16 


Oct. 7, 1969 

6.36 

17.64 

Fall 

Nov. 25, 1969 

7.46 

16.54 

Winter 

Dec, 23, 1969 

7.69 

16.31 

Winter 

Feb. 5, 1970 

7.09 

16.91 

Winter 

Mar. 17, 1970 

6.09 

17.91 

Spring 

Aprxl 28, 1970 

5.04 

18.96 

Summer 

June 9, 1970 

4.36 

19.64 

Summer 

Sept. 1, 1970 

5.45 

18.55 

Fall 


A 

A 


-4 . 7x10 , 
5.8x10 

1 . IxlO 

3.00 

-4 . 4xlot 
5.3x10^ 

0.9x10^ 

1.40 

4 

4 


-3 . 7x10^ 
6,6x10^ 

2 . 9x10 

0.28 

-2 , 1x10^ 
4,4x10^ 

2.3x10^ 

0.23 

-4.5x10^ 

11.7x10 

7.2x10^ 

0.13 

-4. 7x10 J 
6.0x10 

1.3x10^ 

1.60 

-4.8x10^ 
6 . 7x10 

1.9x10^ 

1.00 

-4.2x10^ 

6.5x10 

2.3x10'^ 

0.53 
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were in each case obtained by extrapolation. This "undesirable 
feature, combined with the coarse height and local time resolution 
of the St, Santin facility, combine to make the St. Santin results 
somewhat uncertain, as are those obtained from the Arecibo 
measurements . 

Since the boundary values are important in controlling both 
the magnitude and direction of the terminator fluxes, it is 
worthwhile to examine whether the reversal in direction of the net 
flux that occurs between the Arecibo and St, Santin results can be 
anticipated by comparison of the boundary densities and tempera- 
ture parameters in Tables A.l and A«2« An examination of the 
tables reveals a positive correlation between the 120 km densities, 
when corrected for variations in the temperature parameters, and 
the corresponding computed terminator fluxes. Also, while the 
computed exospheric temperatures in Table A.l are consistently 
higher near sunset than close to s"unrise, the reverse is true for 
the 120 km temperature and it is these latter values, along with 
the boundary densities, that largely determine the magnitude of 
the height integrated energy fluxes. These same trends in both the 
exospheric and 120 km temperatures at Arecibo have been observed 
by Swartz and Nisbet (1971) using data from one summer day 
(June 26) and two winter days (Jan. 13, 14) in 1968. A gradual 
decline in the 120 km temperature from the morning to the evening 
hours is readily apparent from the summer data and is suggested 
by the winter results also, although not as conclusively. 

Fontanari (1973) has performed a harmonic analysis of data 
from St. Santin obtained during equinox and solstice periods for 



63 


the years 1969 and 1970. His results show the height variation of 
the mean temperature for each season along with the amplitude and 
phase of the diurnal temperature variations and reveal enhanced 
temperatures at the evening terminator during summer solstice and 
at both equinoxes. These findings tend to confirm the 120 km 
temperatures summarized in Table A. 2 and together with the results 
of Swartz and Nisbet (1971), point up the latitudinal dependence of 
the lower boundary temperature at the terminators. 

In arriving at global flux estimates under equinoctial condi- 
tions, the geometry of the problem allows the use of symmetry in 
the latitudinal dependence of the energy flux between the northern 
and southern hemispheres. Using the equinox flux values appropriate 
for the latitudes of Arecibo and St. Santin, this assumption provides 
four points which represent the variation of the energy flux with 
latitude. A similar scheme is employed using the solstice results 
by using the values appropriate for the latitudes of St. Santin and 
Arecibo in the northern hemisphere and applying these results to 
the conjugate hemisphere six months later. 

Using this approach and the maximum and minimum values 
quoted previously, the energy flow Q_ into the nighttime thermo- 

1j 

sphere via the low and mid latitude thermospheric winds was found 
to lie between -3.9X10^^ W and 11.2X10^^ W at equinox and between 
-1.5X10^° W and 9.4X10^^ W at solstice. 

3 . 2 The Polar Thermospheric Wind System 

The complexity of the polar thermospheric wind system has 
been deduced from observations and has been considered in detail 
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by a number of authors (Fedder and Banks,, 1972; Meriwether et aL , 
1973; Nagy et al. , 1974; Brekke and Doupnik, 1973; Maeda, 1976). 

The actual wind system is thought to be a superposition of the 
solar tidal winds and the winds due to the motion of the ions which 
are, in turn, influenced by a general dawn to dusk electric field. 
The resultant pattern of the horizontal wind field changes with 
height, since it depends on the degree of coupling between the ion 
and neutral motions . This interaction between the ions and neutrals 
allows a transfer of the electric field energy (of magnetospheric 
origin) to the neutral atmosphere with a resulting temperature 
enhancement and an overall outflow of gas from the polar region 
at heights above 120 km (Brekke et al. , 1973; Nisbet and Glenar, 
1977). 


3.2.1 Method 

The polar region is treated here as an independent energy 
source which supplies a flow of the neutral gas at elevated 
temperatures to the low and mid latitude nighttime thermosphere 
via the south component of the horizontal polar wind system. The 
global energy fl\ix transported by the neutral wind from the high 
latitude regions into the mid latitude thermosphere can be obtained 
using equations (2.13) and (2.14) where in this case the integral is 
over the surface separating the polar caps and the nighttime 
thermosphere, as shown in Figure 2, 

Letting dt be an increment in local time, and assuming 
symmetry in the wind system between the north and south polar 
regions, we may write the energy- input to the nighttime sector as 
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6 

(2)(it/ 12) t) sinedt 

18 



where 


v->- 1 ip 

i 120 ' 


+ 1 n.kT + n. m. g(h-120) dh 


J 


',-(3.7) 


In the above relations Vj^ is the meridional component of the 
horizontal neutral wind directed into the nighttime region. In 
addition, the variation of R with height h is neglected and is 
taken to be the radius at the lower boTxndary (6490 km). The first 
term in equation (3.9) represents the outflow from the polar ■ 
regions between 18 and 6 hours local time. The second term 
represents the energy leaving the system at the lower (120 km) 
boundary. Here, is the magnitude of the mass flow leaving the 
nighttime system and by the mass continuity equation (2.16), the 
defining equation for this quantity can be expressed as 



6 

18 


where Pjyj;(9’t) 



i 120 


(n. m. V, ,) dh 
' 1 1 M 


(3.8) 


( 3 . 9 :) 


The expressions ^) and Pj^(9»t) represent the meridional 

component of the energy and mass flow entering the system at 
geomagnetic cola'titude ±9, and local time t. Similar expressions 
may be obtained for the zonal components using the corresponding 
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zonal component of the horizontal wind. 

An estimate of these quantities, then, requires both the 
direction and magnitude of the neutral wind velocity vector as a 
function of height above 120 km and this has been established using 
the results of 26 high latitude vapor trail measurements. For the 
purposes of this paper, the pattern of outflow is examined using 
geomagnetic latitude and geomagnetic time as the coordinate 
variables. Table 4 shows the geomagnetic latitudes of each of the 
flights, ranging from about 59°N to 74°N with a mean latitude of 
about 67°N. Since all of the information is contained in this 
latitude region, the boundary defining the emitting region has been 
taken to lie at 70°N latitude (0= 20°). 

3.2.2 Boundary Conditions 

The results of each flight provide information on the height 
distribution of the horizontal neutral wind at a fixed local time and 
fixed geomagnetic north latitude near 70°. Using suitably chosen 
boimdary conditions, then, it is possible to estimate the meridional 
components of the energy flow 6j^(9,, t) and mass flow 
each of the flights using equations (3. 7)' and (3'.9).'. 

Rees (1971) has analysed fluorescent emissions from aluminum 
oxide (AlO) and barium oxide (BaO) during a number of twilight 
rocket launches at high latitudes. The results indicate that a 
correlation exists between the morning and evening twilight 
temperatures at 140 and 165 km, and the C9 daily magnetic 
character. This relation is linear in the magnetic index and is 
expressed quantitatively by Rees as 



Table 4. High Latitude Vapor Trail Measurements 


No. 

Location 
N.Lat. E.Long. 

Time 

UT 

Date 

Day Year 

Geomagnetic 
GLT ■ G ^ Lat « 

T 

K 

P 

V 

10.7 

Altitude 
Range (km) 

Ref. . 

1 

65.1 

212.5 

16.7 hr 

288 

72 

4.6 

71.2 

3.0 

95 

110-160 

A 

2 

65.1 

212.5 

15.2 hr 

288 

72 

3.1 

71.2 

3.0 

95 

100-180 

A 

3 

65.1 

212.5 

22.1 hr 

23 

76 

9.6 

71.2 

3.4 

72 

105-125 

B 

4 

65.1 

212,5 

9.6 hr 

86 

73 

21.6 

71.2 

0.2 

74 

95-160 

C 

5 

65.1 

212.5 

7.7 hr 

50 

74 

20.0 

71.2 

.5 

85 

95-132 

C 

6 

58.8 

265.7 

3.5 hr 

142 

63 

21.8 

73.4 

.7 

91 

80-175 

D 

7 

58.8 

265.7 

4.1 hr 

143 

63 

22.4 

73.4 

1.2 

89 

80-160 

D 

8 

58.8 

265.7 

8.3 hr 

143 

63 

2.2 

73.4 

.6 

89 

85-160 

D 

9 

58.8 

265.7 

6.0 hr 

306 

64 

0.4 

73.4 

1.2 

75 

95-130 

D 

10 

58.8 

265.7 

0,0 hr 

58 

65 

18.0 

73.4 

1.8 

74 

90-160 

D 

11 

67.8 

20.4 

7 . 7 hr 

32 

68 

10.9 

58.5 

1.6 

262 

90-160 

E 

12 

67.8 

20.4 

7.6 hr 

35 

68 

10.8 

58.5 

3.8 

111 

90-165 

E 

13 

67.8 

20.4 

18.5 hr 

289 

68 

21.3 

58,5 

1.3 

138 

90-155 

E 

14 

67.8 

20.4 

6.7 hr 

306 

68 

9.8 

58.5 

4.3 

156 

90-170 

E 

15 

80.6 

58.0 

13.6 hr 

57 

70 

18.3 

65.9 

2.6 

182 

100-160 

F 

16 

80.6 

58.0 

16.9 hr 

76 

71 

21.7 

65.9 

2.2 

117 

100-160 

F 

17 

80.6 

58.0 

17.2 hr 

77 

71 

22.0 

65.9 

1.1 

115 

100-160 

F 

18 

80.6 

58.0 

10.2 hr 

307 

• 71 

15.2 

65.9 

.8 

117 

100-160 

F 

19 

80.6 

58.0 

9.3 hr 

311 

71 

14.4 

65,9 

2.2 

107 

100-160 

F 

20 

80.6 

58.0' 

12.9 hr 

52 

72 

17.7 

65.9 

2.4 

194 

100-160 

F 

21 

80.6 

58.0 

16.9 hr 

76 

72 

21.7 

65.9 

2.7 

135 

100-160 

F 

22 

80.6 

58.0 

21.8 hr 

79 

72 

2.4 

65.9 

1.2 

136 

100-155 

F 

23 

80.6 

58.0 

12.3 hr 

46 

74 

17.1 

65.9 

1.1 

86 

100-140 

F 

24 

80.6 

58.0 

17.5 hr 

78 

74 

22.3 

65.9 

.2 

81 

100-140 

F 

25 

80.6 

58.0 

16.3 hr 

72 

71 

21.0 

65.9 

4.9 

114 

110-160 

F* 

26 

80.6 

58.0 

22.9 hr 

72 

71 

3.5 

65.9 

5.0 

114 

115-160 

F* 

References : A. 

Bedinger 

(1973) , ^ 

Nielsen 

(1976), 

G. Bedinger (1975) 

> Ri. 

Bedinger 

(1966), 



E. Rees (1971), F. Andreeva et al. (1975). 


tKp values represent average of 12 hrs. preceding flight. 
*Data obtained from six spherical clouds. 



T(140 km) 


420 + 25 • C9 


K 


(3. IT)) 


T(l65km) = 450 + 50 • C9 K {3.11) 

for morning twilight and 


T(140 km) = 

500 + 25 • 

C9 

K 

(3.12) 

T(165 km) = 

550 + 50 • 

C9 

K 

(3.13) 


for evening twilight. 

In the absence of any additional experimental information, it 
is assumed here that the temperature change between the evening 
and morning sectors at these altitudes varies linearly with local 
time in the nighttime hemisphere. This assumption allows the 
temperatures at 140 and 165 km to be estimated at each of the 
launch times for the rocket firings previously discussed. For 
purposes of this study, then, the adopted temperatures between 18 
and 6 hours become 


140 ^ 

540 - 6. 67 mod (t 4- 12, 24) 

+ 25 • C9 

K 

(3.14) 

165 

600 - 8.33 mod (t + 12, 24) 

+ 50 * C9 

K 

(3.15) 


where t is the geomagnetic local time in hours. 


Additional boundary parameters to be employed are the 
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120 km atomic oxygen and molecular nitrogen densities and the 
exospheric temperatures at the times and locations of each of the 
rocket flights listed in Table 3. The OGO-6 satellite quadrupole 
mass spectrometer has been used to measure O and N 2 densities 
in the height range from 400 to 600 km (Hedin et al, , 1974). In 
addition, neutral temperatures in the altitude range from 200 to 320 
km have been obtained using 630 nm oxygen emission line results 
from the Fabry-Perot interferometer on board the OGO-6 (Blamont 
and Luton, 1972). The neutral temperatures at these altitudes are 
within 100 K of the exospheric temperature and a good approxi- 
mation of T^ is obtained by using the method of Nisbet et al. (1973). 

The vertical temperature distribution above 120 km can then 
be found using this value for T^ and the temperatures at 140 and 
165 km obtained from equations (3.17) and (3.18). Using these ' 
values, equation (3.4) can be solved to obtain 


s 





(3'. 16) 


and 


120 


= - (T^- T^4q) 


^140 

"^03 ■ 


-1 4/5 


(3.17) 


The results of the OGO-6 density and temperature measure- 
ments have been grouped according to latitude, season, geomag- 
netic local time and Kp range. Measurements made between 65° 
and 75°latitude and near spring and fall equinox are particularly 
useful, since the polar boiindary is fixed at 70°N latitude and since 
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the symmetry of the problem during equinox allows the use of data 
from both hemispheres. The local time variation of the densities 
and temperatures is examined by collecting and averaging all data 
obtained in increments of 1 hour. The time variation can then be 
represented by a series of 24 values, from 0 to 24 hours 
geomagnetic local time. 

In using these results to estimate the atomic oxygen densities 
at 120 km, the diffusive equilibrium assumption is first made and 
the effects of departures from diffusive equilibrium are then taken 
into account, Nisbet and Glenar (1977) have shown that the atomic 
oxygen densities at 120 km and near 70^N geomagnetic latitude can 
be related to the diffusive equilibrium values by 

logjt, n(O)i 20 = logj(| n( 0 )* 2 j + 0. 96 (3.18) 

where n(O)j ^20 estimate of the 120 km atomic oxygen density 

taking into account the effects of diffusion 
n(O )^20 diffusive equilibrium value 

v^( 0 )^ 2 q is the average upward vertical velocity of atomic 

' 

oxygen at 120 km, in msec 

Using horizontal winds obtained from the vapor trail 
measurements listed in Table 4, Nisbet and Glenar, have derived a 
value for ’'^2^^^20 0,76 msec appropriate for geomagnetic 

latitudes close to 70 degrees. Inserting this value in equation (3,I8)„ 
the required correction for the atomic oxygen densities at 
geomagnetic latitudes near 70 degrees becomes 
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n(O)^20 = 

Figure 6 shows the boxmdary conditions calculated in the 

I 

above manner using the temperature results of Rees (1971) and the 
OGO-6 temperature and density data. Since the rocket firings 
occur xinder a variety of geomagnetic activity levels, boundary 
/parameters have been determined for three ranges of Kp: 0,0-1, 2, 
1.3 -4,0 and 4,1- 9.0. In each case one sees an enhancement of 
both the temperatures and the atomic oxygen densities near the 
evening terminator and a corresponding depression near the 
morning terminator. 

3.2.3 Results of the High Latitude Analysis 

Analyses of the motion of rocket vapor trails released at high 
latitudes show the height variation of wind velocity and direction 
at the location of each firing. Using these results, and values of 
the boundary temperature and density profiles appropriate to each 
flight, equations (3.6) and (3.8) have been used to estimate the 
height integrated horizontal energy and mass flux above 120 km in 
each case and the resulting net energy exchange from the high 
latitude region to the nighttime thermosphere. ' The results from 
three of the twenty six flights listed in Table 4 have not been used 
in evaluating these fluxes. Wind velocity data for flight number 3 
was only obtained to an altitude of 125 km. Also, the results of 
measurements 25 and 26 were obtained from observations of 
spherical clouds during magnetically disturbed periods with a Kp 
index at flight time of approximately 5. The choice was made to 



Figure 6a Variation with local time of boundary parameters during spring and fall equinox for 
Kp 0.0 - 1.2. Data is representative of 70°N and 70°S geomagnetic latitudes. 




Figure 6b Variation with local time of boundary parameters during spring and fall equinox for 
Kp 1.3 - 4.0, Data is representative of 70°N and 70°S geomagnetic latitudes. 


ozt 




Figure 6c Variation with local time of boundary parameters during spring equinox for 

Kp 4.1 - 9.0. Data is representative of 70 H and 70°S geomagnetic latitudes. 
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restrict our analysis to data representative of only relatively quiet 
magnetic conditions (Kp<4.0), Tbe results of these measurements, 
then, were also excluded. 

Figures 7 and 8 show the height integrated energy and mass 
fluxes above 120 km as determined for each of the 23 high latitude 
rocket flights. As is expected, the flow patterns in both cases are 
similar, while slight differences in the relative lengths of the 
arrows between the figures are the result of weighting effects from 
the different height distributions of the parameters involved. 

Most of the flights observed here were in or near the 
statistically determined auroral belt (Heppner,, 1372)^ Accordingly, 
many of the features seen in Figures 7 and 8 have been observed 
previously using results from a number of high latitude barium 
cloud releases (Meriwether et al., 1973) and simultaneous 630 nm 
airglow and incoherent scatter observations (Nagy et al. , 1^74). 
Figures 7 and 8 show; 

(i) Large, southerly directed flows in the predawn sectbr.. 

(ii) Reversal of the zonal component near 2300 hours 
geomagnetic local time. A shift occurs here from a 
westward zonal component in the premidnight sector to a 

j 

generally eastward drift in the morning sector. 

(iii) An anomalous west to east drift for the flights in the 
prenoon sector, where the effects of the diurnal solar 
tide would be expected to have the greatest influence and 
so give rise to an east to west drift component away from 
the afternoon diurnal bulge. The three flights used to 
provide flvix estimates in the prenoon sector have been 
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SCALE 

t * *-t « M ^ 

O 3 6x10 ‘ 

kg irf^stcT* 


18 hrs 



Ohrs 


GEOMAGNETIC LOCAL TIME 


Figure 7 


I 


Total mass flow above 120 km derived from 23 high latitude 
vapor trail measurements . The center of each arrow 
corresponds to the location of each measurement. 
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SCALE 

I ■ 1 ' I ^ 

0 1 2 3x10^ 

Wm"' 


18 hrs 



0 hrs 

GEOMAGNETIC LOCAL TIME 


Figure 8 Energy flow above 120 km derived from 23 high latitude 
vapor trail measurements. The center of each arrow 
corresponds to the location of each measurement. 
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taken from Rees (1971), who has found a correspondence 
between the west to east neutral winds and south directed 
ground magnetic fields, which are present during each of 
the prenoon rocket flights employed here. No relationship 
of any kind has been found relating the ground magne- 
tometer data and the north- south neutral wind, however. 
Thus, it has been assumed that the north components 
obtained here are representative of that geomagnetic 
latitude and local time. 

A correlation has been found to exist between the local time 
dependence of the meridional components €^^(0", t), Pj^(6, t) and that 
of the 120 km atomic oxygen boundary densities, plotted on a 
logarithmic scale (Nisbet and Glenar, 1977), These relationships 
are shown graphically in Figures 9 and 10, Data for Kp 0,0-1, 2 
has been used here, although the same correlation exists for 
results from the other Kp ranges as well. Since data on the 
boundary densities are available over the entire 24 hours, this 
relationship allows one to estimate the meridional flows between 
5.0 and 12.0 hours geomagnetic local time where there is a 
paucity of experimental wind results. Linear regression analyses 
between the densities and the fluxes were performed and the 
following relationships found. 

= 1.638X10^ log^Qn(0)^2Q - 2.806X10^ W m~^ 

(3,20) 

kg m"^sec"^ 


(3.21) 



TOTAL MASS FLUX (kg tTr*se<f*) 


Figure 9 



North components (g^(0,t)) of total mass flux at 70°N geomagnetic latitude (0 = ! 
height integrated above 120 km. values of the derived 120 km atomic oxygen 

density (Kp 0.0 - 1.2) have been shown for comparative purposes. 


log,pn(0) 




COMPUTED ENERGY FLUX 


17.20 


24 2 4 

DCALTJME (hr) 


d K) 


1 energy flux at 70°N geomagnetic latitude (5 
values of the derived 120 km atomic oxygen 
own for comparative purposes. 
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The primes in equations (3.20) and (3.21) have been included to 
distinguish between the experimental values and the best fit values 
obtained from the regression analysis. Positive values of these 
quantities are those directed to the north in the northern hemi- 
sphere. Using these relations, equations (3.6) and (3. 8) have been 
solved to obtain a value of approximately 2.3X10^^ W (23 GW) as 
the energy input from the north polar region into the system at 
70°N geomagnetic latitude (0= 70°). Since symmetry has been 
assumed between the high north latitude and the conjugate south 
polar regions, the estimate arrived at for the net energy input 
Qp to the nighttime system is, for Kp values from 0.2 -4.0, 

Qp = 4.6 ± 2.3 X 10^° W (3.22) 



CHAPTER IV 


ENERGY LOSS PROCESSES 

4.1 Downward Conduction at the Lower Boundary 

The major heat loss from the night hemisphere is downward 
conduction due to the large vertical temperature gradient. The 
important parameters in establishing absolute magnitudes of the 
conduction loss are the temperature gradients, and to a lesser 
extent the absolute value of the temperature and the relative 
densities of the atomic and molecular major constituents which 
affect the conduction. 

From equation (2»15), the downward conduction at the 120 km 
boundary is given by 

T 8T T -2 -1 , 

C = A J XXI sec (4.1) 

Here, 9T/9z is the vertical temperature gradient, and X. is the 
weighted mean thermal conductivity for the gas mixture. The values 
used for the thermal conductivity of N^, and O are those given 
by Bauer (1972) and are expressed as follows 


X(02) = 

1.86 X 10"^ 

^ .84 
■‘‘n 

(4.2) 

X(N2) = 

2.72 X 10“^ 

T .80 

n 

(4.3) 

II 

O 

6.71X 10'^ 

.j, . 71 
■^n 

(4.4) 


The mean thermal conductivity can then be expressed as 
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X 


XjN^) + X(02) n(02)/n(N^) + X(0) n(0)/n(N2) 
1 + n(02)/n(N2) + n(0)/n(N2) 


(4.5) 


Using the 1972 Cospar International Reference Atmosphere (CIRA, 
1972), representative density ratios at 120 km have been evaluated 
and used to compute the mean thermal conductivity. The result is 

X = 1.8X10 T • + 1. 72X10 T • !*1- » 

n n 

+ 1.64X10“^ T * J m"^ sec"^K"^ (4.6) 

n 

A number of difficulties are encountered in making experi- 
mental measurements of thermospheric temperatures at night. 
Probably the best estimates of temperature in the region above 
120 km are obtained from incoherent scatter (Wand, 1974; Swartz 
and Nisbet, 1971). Unfortunately, electron densities at night in the 
region are so low that up until now only daytime measurements at 
these heights have been available. Individual rocket measurements 
of temperature in this region show a pronounced wave structure 
with considerable day to day variation with amplitudes of the order 
of 50 K and a phase variation with height. These factors combine 
to make individual in- situ measurements rather difficult to relate 
to average nighttime energy fluxes. 

While the incoherent scatter data do not provide nighttime 
temperature gradients, it appears that the dominant tidal component 
is semidiurnal rather than diurnal (Salah and Wand, 1974) so that 
the uncertainty should not be large. In accord with the analysis by 
Fontanari (1973) of the temperature amplitudes and phase 
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distributions and their associated errors (Section 3.1.4), the height 
and local time dependence of the temperature at St, Santin is 
given by 

T(t, z) = T (z) + A(z) cos {-jr/12)(t - <#),(z)) 

O Cl 

+ B(z) cos (ir/6)(t- ^ (z)) (4.7) 

where T^(z) isl- the seasonal mean temperature 

A(z) is the amplitude of the diurnal component 

B(z) is the amplitude of the semidiurnal component 

^^(z) is the phase of the diurnal component in hours 

<l> (z) is the phase of the semidiurnal component in hours 

s 

The average value of the conduction over the night hemisphere 
at the latitude of St. Santin can then be obtained using equations 
(4,1), (4.6) and (4.7) and is given by 

t 

r 

= - 24^ (/-t ) y \9T/az dt (4.8) 

® ^ t 

s 

where t and t are the times of local sunrise and sunset at 
r s 

St. Santin for the season under consideration 

The conduction rate associated with the steady or diurnal average 

temperature gradient was calculated to be 5.2X10"^ Jm~^'sec"^ 

-4 -2 -1 

for spring equinox, 4. IX ID Jm sec for fall equinox and ' 

- 4 - 2 -1 

4,8X10 J m sec for summer solstice. The corresponding 



85 


ranges of solar Fj^q f lux were 123-232, 102- 207 and 109- 239 
-2 

W m. Hz for spring equinox, fall equinox and summer solstice 
respectively. No analysis was performed for winter solstice, since 
the results of the measurements made were inconclusive. While 
there is considerable uncertainty associated with the diurnal 
component, the 120 km conduction averaged over the nighttime 
hours is found to be lower than the seasonal mean in each case 

with considerable variation in the amount of reduction. The values 

-4 

of average nighttime conduction were found to be 3.7±1.4X 10 , 

-4 -4 -2 -1 

3. 4 ±1.6X10 and 4.7±.7X10 Jm sec for spring equinox, 

fall equinox and summer solstice respectively. 

Similar analyses have been made at Millstone (42.6°N) by 

Salah (1974) and at Millstone and St. Santin (44.6°N) by Salah and 

-22 -2 -1 

Wand (1974) for values of F^q ^ between 95 and 186 X 10 W m Hz . 

These studies provide estimates for the diurnal average conduction 

_4 _2 -1 

of 5.7 and 6.6 X 10' Jm sec respectively. 

Based on a statistical analysis of measured temperatures 
between the years 1966 and 1970, Rohrbaugh,(1972) has analysed the 
local time variation of the 120 km conduction at Arecibo. The 
results are represented by the Fourier series 



(4.9) 


where the sets of coefficients A. and B. have been determined by 

J J 

Rohrbaugh for equinox and winter and summer solstice conditions. 
When the steady state temperature given by the Aq term is used to 
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determine the thermal conductivity, the average nighttime conduc- 

_4 

tion calculated using equation (4.8) is found to be 2. 8 ±.8X10 

-2 -1 . -4 -2 -1 

Jm sec for spring and fall equinox, 2.3±.5X 10 Jm sec 

-4 -2 -1 

for summer solstice and 2.2 ±,9X10 Jm sec for winter 

-22 

solstice, for values of F^q y flux between 100 and 225 X 10 

W m“ Hz” . The xmcertainty in each of the above values has been 

deduced from statistical uncertainties as well as from the errors 

in the temperatures near 120 km, as deduced from the Arecibo data. 

Figure 11 shows the results derived from both the St, Santin 

and the Arecibo data and indicates a marked difference in both the 

steady state conduction and the amotint of variability between the 

two stations. For comparative purposes, individual values of 

conduction derived from isolated pitot tube measurements at 

Arecibo, Wallops Island and Ft. Churchill have been included. 

These data were obtained for values of Fj^q y between 108 and 174 
-2 -1 

W m Hz . A discussion of the method of obtaining temperature 
gradients from pitot tube measurements has been given by Smith 
et al. (1967). 

The variability of these individual measurements is assumed 
to be the effect of reversible heating from atmospheric gravity 
waves which exhibit a characteristic phase variation with height. 
Hines (1965) has shown that the amplitude of temperature fluctua- 
tions increases from roughly 10 K low in the E-region to as much 
as 30 K near 110 km. In addition, from an analysis of Na cloud 
drifts, Kochanski (1964) has deduced vertical wavelengths of the 
order of 16 km near 120 km, thought to result primarily from the 
effects of gravity waves .. Assuming a temperature amplitude of 
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Figure 11 Diurnal variation of downward heat conduction at 120 km. 
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30 K and a representative vertical wavelength of 18 km, a maximum 

departure from the steady state temperature gradient of about 

10 K/km is found to occur. At a 120 km temperature of 350 K, 

-4 -2 -1 

this corresponds to deviations of the order of 3X10 Jm sec 
from the steady state conduction at 120 km. The size of this 
deviation is represented by the error bars in Figure 11, 

If the latitudinal dependence of the nighttime conduction rate 
is adequately known, the total conduction over the 120 km boundary 
of the nighttime thermosphere for a given season is expressed as 

70° 

= 1.85X10^^ \ <C> L, cos<p d<P W (4.10) 

-70° 

where L is the length of night in hours 

- z 

C > is the average nighttime conduction rate in W m” 

<p is the latitude in degrees 

As in the treatment of~the horizontal winds, the geometry at 
equinox allows the assumption of latitudinal symmetry in the 
conduction rate between the northern and southern hemispheres. 
Using the equinoctial conduction rates at the latitudes of Arecibo 
and St. Santin, four points are available which represent the 
variation of the conduction with latitude. This allows an approxi- 
mate solution of equation (4.10) and yields a value of 

= 8.6 ± 2, 7 X 10^^ W equinox (4.11) 

for the total downward conduction at times close to equinox. 
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For either summer or winter solstice conditions, the latitude 
dependence of the conduction is again obtained by using the rates 
appropriate for the latitudes of St. Santin and Arecibo in the 
northern hemisphere and applying these results to the conjugate 
hemisphere six months later. 

Since it has not been found possible to estimate the winter 
solstice conduction rate at St. Santin, however, two cases have 
been studied. The first case assumes that the conduction rate is 
the same at winter solstice as it is at summer solstice. The 
second assumes that the rate oscillates aroxuid the equinox value 
with a 12 month period so that based on the results obtained, the 
winter solstice value lies below either of the equinox values. When 
this approach is taken, the following values for the total conduction 
result. 


11 

o 

a 

8.1 ± 2.6 

X10^° 

W 

solstice 

maximxim 

(4.12) 

11 

o 

a 

6.6 ± 2.1 

X10^° 

W 

solstice 

minimtim 

'(4.13) 

Infrared Radiative 

Losses 





4.2.1 

Emission at 

63 u from 

Thermospheric Atomic 

Oxygen 


Kockarts and Peetermans (1970) have discussed the contri- 
bution to thermospheric cooling by infrared radiation from atomic 

oxygen in the height range between 50 and 250 km. The ground 

3 3 

state of atomic oxygen consists of three low lying levels Pq, P^ 
3 

and kinetic energy of motion of the constituents is 

converted into infrared radiation by collisional excitation and 

3 3 

subsequent spontaneous deexcitation to either the P^ or the P^ 
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3 3 

levels. As shown in Figure IZ, the levels and lie .028 and 

3 

,020 eV above the ground level P^ and emit at 147 fx and 63 jx 
respectively. The transition probabilities A., are simply the 
spontaneous transition rates per atom and are listed by Wiese et 
al. (1966) as A^^= 1. 70X 10“^ sec"^ and A^^ = 10"^ sec^^ The 

energy lost per unit volume from this mechanism may then be 
written 


^ " ^63(x ^147 PL 

^ “ ^1-^12^1^^^ ^^0 " ®1^ "^Ol 


«. 14a) 
(4. 14b) 


3 3 

where and are the energies of levels P^ and P^^ with 

3 

respect to P^ (here, ^2 taken to be zero). 

Ay^ and A^^^ are the spontaneous emission coefficients from 
levels 1 to 2 and 0 to 1 respectively 
nj^(O) and iiQ(O) are the populations of levels 1 and 0 
respectively 


In the region of interest here, the collision frequency between 
atomic oxygen and all constituents is ass\omed high enough to 
maintain a Boltzmann distribution, which relates the population of 
level i, i.e, n.(0), to the total number of atomic oxygen atoms. 

This relation is given by 

g exp(-e,/kT) 

n (O) = n(0) (4.15) 

S g. exp (-£ ./kT) 

j ^ J 
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where g. and are the statistical weights of levels i and j 
respectively 
T is the temperature 

n( 0 ) is the total atomic oxygen density 


At temperatures of the order of 500 K, nearly all atoms are in 

one of the above three levels. Using values for the statistical 

weights of gQ=l> ^ volume emissions ^ 1^3 

and L, are foxmd to be 
147 |i 


1.69X10'^^ n(0) exp'(-228/T) -Wm“^ 

1 + 0,6 exp(-228/T) + 0.2 exp(-326/T) (4.16) 

4.59X10"^'^ n(0) exp(-326/T) Wm"^ 

1 + 0 , 6 exp (-228/T) + 0. 2 exp (-326/T) (4.17) 


where n( 0 ) is the number of oxygen atoms per m 


-3 


Inspection of the above equations quickly shows that for the 
same values of T and n(0), the 147 [x contribution to atomic oxygen 
emission is negligible. In consequence, only emission at 63 jx will 
be considered. 

The emission rate given by equation (4.14) applies only for an 

optically thin atmosphere, where the radiation travels unimpeded 

and ultimately leaves the atmosphere. In reality, some absorption 

of the radiation occurs in the reexcitation of atomic oxygen from 

3 

the level and the rate of this process is directly proportional 

3 

to the population of the level. As a result, the emission at a 
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given height is reduced and the departure from optically thin 
emission increases as height decreases. Kockarts and Peetermans 
have obtained a reduction factor y to be applied to equation (4.14) 
expressed as a fraction of the emission. Their results are shown 
in Figure 13 where curves corresponding to exospheric tempera- 
tures of 750 K and 2000 K have been included. The average of 
the two curves shown was chosen to represent a 1200 K exospheric 
temperature and used to compute the height profile of the emission. 
The resulting curve was fitted using the least squares technique 
and a good representation for z^l20 km is given by 

1-X = 1200 K (4.18) 

where z is the height in km. 

For the purpose of estimating radiative losses, a first order 
approximation to the atomic oxygen density distribution above 120 
km is obtained assuming diffusive equilibrium. Using the tempera- 
ture distribution given by equation (3.4), the height distribution of 
atomic oxygen (z in km) is thus given by 

n^(0) = T^ 

(4.19) 

with z = 120 km. 

o 

The emission (q/ « ) from a column of unit cross section and 

'^63 (jl' 

extending upward from 120 km can then be calculated using 


iriog 

kT„ 


z-z^+(l/s)ln(T^/T^2,) 



ALTITUDE (km) 
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Figure 13 Reduction factor resulting from radiative transfer applied 
to 63p emission. (Reference: Banks and Kockarts, 1973) 
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120 

where ^“X determined using equations (4,14) 

through (4.19) with appropriate boundary values. 

Donahue et al. (1973) have determined the distribution of 

atomic oxygen as a function of altitude and latitude from OGO - 6 

observations of the atomic: oxygen night air glow emission. 

They find that at 120 km, a vadue fojr the ^atomic .oxygen -de'nsity 
17 -3 

of 1.1X10 m represents all but the most unusual conditions, 
while the observed values ranged from 0.41X10^^ to 1. 35 X 10^^ m"^ 
(Donahue, 1976). It seemed appropriate, then, to assume that the 

120 km atomic oxygen density lies usually between 0.9 and 

17 -3 

1.2X10'-^ m 

In order to establish a representative temperature profile, 
temperature parameters derived from the September 17, 1974 
Arecibo incoherent scatter data (Table A.l) have been employed. 
This day represents low geomagnetic activity and the scatter 
results obtained on this date were among the best of the available 
measurements. The extremes of temperature listed in Table A.l 
for this day should represent the maximum and minimum tempera- 
ture over the entire 24 hours since the temperature oscillations 
at these heights are predominantly semidiurnal and since data was 
collected over a period of about 12 hours. 

Using the above extremes in the density and temperatures, 
equation (4.20) has been evaluated and the 63 p. loss rate per square 
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meter was found to be 


^63 [X 


0.86 ±0.13 X 10 


-4 


W m 


-2 


(4,21) 


4.2.2 Emission at 15 |jl from Thermospheric Carbon Dioxide 
An additional source of infrared radiative loss above 80 km 
arises from radiative deexcitations of the bending vibration band 
of CO^ at 15 . The basic considerations involved in calculating 

the radiative loss rate from this process -are given by Houghton 

- 6 

(1969). For total pressures significantly less than 10 atmospheres, 
collisional deexcitation of the excited states in the band becomes 
negligible compared to radiative deexcitation and in addition, the 
medium becomes optically thin for this transition. The radiative 
loss rate per unit volume over the entire band is then given by 



hy- So 

n(CO.,) — 

2' T g. 



-1 


(4..2'2) 


In equation (4.22) hv is the average transition energy over the 

-13 

band (hv'^l, 33X10 erg) while and g^ are the statistical 

weights with values of 1 and 2 respectively. .The estimate obtained 

by Houghton for the vibrational relaxation time t at 210 K is 

c 

6.0X10 ^/p sec where p is the total pressure in atmospheres. 

This value will be used in this analysis, although Houghton 
indicates that the value of decreases with increasing temperature. 


The term (g^/gj^) 


exp (hv /k T ) - 1 


"1 -1 


is the ratio of the probabilities 


of excitation from the groimd state to deexcitation from an excited 
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state when allowance is made for all excited states that comprise 

the band. Assuming a height independent volume mixing ratio 

-4 

for CO^ of 3.3X10 (Houghton,, 1.9'69,), equation (4.9) becomes 


15 p 


_4A 

1.99X10 


1 2 


n{W^) + n{0^) + n(0) 



W m”^ 


(4% 23)' 


The emission above 120 km from a column of unit area is then 
given by 

OD 

120 


Because of the quadratic dependence on the total density, 
rather large variations in the 15 p emission rate would normally 
be expected. Using the range of 120 km N 2 ^ and densities found 
in Table A,1 for September ,17, along with the previous boundary 
values, the 15 p emission has been found to be 

qi5 p ~ 0.30 ± 0.20 X 10“^ W m"^ (4.25) 

4.2,3 Total Radiative Loss 

From the preceeding sections, the sum of the 63 p and 15 p 
emissions above 120 km from a column of unit cross section is 


^63 p 


+ q 


15 


|J- 


1.16 ± 0.24 X 10 


-4 


W m' 


-2 


(4.26) 
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where the error is determined by the bounds placed on the 120 km 
boundary values, as discussed previously. Over the entire 
nighttime thermosphere, the above value corresponds to a total 
loss rate of 

= 3.0 ±0.6X10^^ W (4.27) 

for the contributions of 63 [jl and 15 jJ- emission by thermospheric 
atomic oxygen and carbon dioxide to the energy balance. 



CHAPTER V 


RESULTS OF THE -ANALYSIS 

5 . 1 Concltisions 

5.1.1 Conductioir and" Radiation Losses 
Two sources of data examined in establishing the lower 
boundary temperature gradients were incoherent scatter and rocket 
pitot tube measurements. From the St. Santin measurements, it 

was found that the mean nighttime conduction flux was 3.5X10~^ 

-2 -4 -2 
W m at equinox and 4.7X10 W m at solstice for values of 

the 10.7 cm solar flux between 102 and 239X10 W m”^ Hz~^. 

The values for Arecibo were measured at times when the 10.7 cm 

-22 -2 -1 

solar flux was between 100 and 225X10 W m Hz and were 
found to be 2.8X10 ^ W m ^ at equinox and 2.3X10~"^W m~^ at 
solstice. 

In the case of the pitot tube measurements it was found 
difficult to make adequate estimates of the average temperature 
gradients because of the large wave components that are always 
present. The pitot tube estimates were found to lie between^3 X 10*"^ 

and 9X10 ^ W m ^ at equinox and between 2X10 ^ and^^lOX lO”^ 

_2 

W m at solstice. 

It appears that the incoherent scatter sounders provide the 
best data for measurements of the daytime temperatures and their 
gradients in the region near 120 km. It is not possible at the 
present time, however, to measure nighttime temperatures in this 
region at mid latitudes using incoherent scatter. Using only 
daytime data to estimate the diurnal component of the temperature 



100 


variation is not at all desirable. The small niimber of incoherent 
scatter stations also makes any assumed latitudinal variation 
conjectured. All that can be said for the pitot tube measurements 
is that they do not show any large diurnal component in either the 
temperatures or the gradients at 120 km and that the values are 
not inconsistent with the estimates based on incoherent scatter 
observations . 

Based on the available data it was estimated that at equinox 
the downward energy conduction over the hemisphere was between 
5.9X10^° and 11.3X10^° W at equinox and between 4.5X10^° and 
10.7X10^^ W at solstice. The radiation loss from the nighttime 
thermosphere from atomic oxygen and carbon dioxide was estimated 
to be between 2.4X10^*^ and 3.6X10^° W. 

5.1.2 High Latitude Thermospheric Winds 

The contribution from the neutral winds blowing in the 
meridional direction was estimated using 23 high latitude rocket 
vapor trail measurements. These were made at different times, 
seasons and latitudes. Because of the complexity of the high 
latitude neutral wind system, a far larger data base would be 
desirable. Still, the general features of the wind pattern appear to 
be understood and the clear correlation with the OGO - 6 density 
data makes it appear that the results are of the right order of 
magnitude. Based on these wind velocities and high latitude density 
data it was estimated that this source provided between 2.3X10^^ 
and 6.9X10^^ W to the nighttime thermosphere. This includes the 
contribution from high latitude heating as well as that from the day 
to night flow across the pole. 
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5.1.3 Low Labitude Thermospheric Winds 

Because of the coordinate system chosen for this analysis 
it is convenient to consider two components separately, one due to 
the winds measured by an observer on the ground and the second 
due to the motion of the terminator with respect to that observer. 
For these components, quite different pieces of experimental data 
control the estimates of the energy. For the winds measured with 
respect to the observer fixed to the earth, the difference in the 
wind profiles at the morning and evening terminators is of prime 
importance. For the component due to the motion of the thermo- 
sphere across the terminator the velocity is known exactly and the 
dominant terms involve the difference in the densities and tempera- 
tures at the two terminators. 

The winds relative to the ground were measured using data 
from 23 pairs of morning and evening rocket flights. Based on 
these measurements alone, the derived energy fluxes ranged from 
0.24X10^ W m"^ to 0.48X10"^'W m“^ of terminator length. 

The preceding term has attracted most attention in estimates 
of energy transport from the day to the night hemisphere. The 
reason for this is that little direct data has been available regarding 
density and temperature variations in the region near 120 km. 

As a result, CIRA 1965 and CIRA 19,72 use approximately constant 
densities, temperatures and temperature gradients at 120 km. When 
changes are incorporated, they reflect primarily attempts to 
compensate for the shortcomings of diffusive equilibrium models in 
explaining the variations in density at satellite altitudes. Because 
of effects due to tidal components and wind induced diffusion, such 
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models cannot be expected to provide good estimates of the relevant 
parameters in the 120 km region. It is thus not possible to use 
them in theoretical investigations of these effects.. 

In the present work it has been shown that the corotating 
winds, rather than the meanured winds, are the dominant contri- 
butor in determining the energy gain or loss to the system from 
the winds, and that there is little density and temperature :data "on 
which to base good estimates of the energy exchange from this 
source. The semidiurnal tidal modes, although they can be well 
determined from incoherent scatter, do not affect the calculations 
in the first order and it is not possible to accurately resolve the 
diurnal amplitude from data obtained only during daylight hours. 

The estimates of the densities near the terminators were found to 
introduce large errors, and while the temperatures at the 
terminators rarely change by more 'than 15 or 20 percent of the 
average value, as observed in Tables A.l and A. 2, the densities 
quite commonly change by a factor of 2 or more. The conclusion 
to be drawn is that a knowledge of the boxmdary densities is of 
even greater importance than the temperatures in establishing the 
terminator energy and mass fluxes. 

5. 2 Summary of Contributions to the Energy Input and Loss 

Processes 

The estimates arrived at for the energy input and loss terms, 
as previously discussed, are listed in summary form in Table 5. 
From a comparison of the inputs and losses, there would appear 
to be a need, under both equinox and solstice conditions, for an 
additional energy source to satisfy the energy balance. The 


Table 5. Summary of Energy Inputs and Losses 


Input Terms 

Average 
(Units of 
10 W) 

Range 

(Units of 10^0 W) 

Loss Terms 

Average 
(Units of 
10^° W) 

Range 

(Units of 10-^^ 

Low and Mid Latitude Winds 




Conduction 



Equinox 

3.7 

-3.9 

11.2 

Equinox 

8.6 

5.9 11.3 

Solstice 

4.0 

-1.5 

9.4 

Solstice 

7.6 

4.5 10.7 

High Latitude Winds 

4.6 

2.3 

6.9 

Radiation 

3.0 

2.4 3.6 


TOTAL; TOTAL; 


Equinox 

Solstice 


8.3 

8.6 


Equinox 11.6 
Solstice 10.6 


o 

OJ 
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average deficit is 3.3X10^^ W (2. OX 10^^ W) for equinox (solstice) 
with a possible range from a 16.5X10^^ W (13.5X10^ W) deficit 
to a 7.8X10^^W (9.4X10^° W) excess. 

It is interesting to compare these estimates of the energy 
deficit with terms that have not been included in the analysis, 

Ching and Chiu (1973) have modeled the volume Joule heating rate 
at local noon under spring equinox conditions and include the day 
to night heating ratio as well. If their rates are considered 
representative of the entire nighttime thermosphere, the total 
Joule heating is found to be 1.7X10^*^ W. 

Using the results of satellite observations of neutral density 

perturbations, Klostermeyer (1973) has selected a dominant gravity 

wave mode and determined the energy dissipation by viscosity and 

ion -drag at thermospheric heights and for geographic latitudes 

near 50°N. By height integrating the calculated heat input above 

-4 - 2 

120 km, a value of 1.7X10 W m is obtained. If this rate is 
applied over the entire nighttime thermosphere, the total heat 
dissipation from this gravity wave mode becomes 4.2X10^^ W, 

These estimates of Joule and gravity wave heating constitute 
a sizable portion of the apparent discrepancy between the energy 
input and loss terms arrived at in this preliminary study. It is not 
proper to considef^he above two sources as completely independent, 
however, since mechanisms have been proposed which couple Joule 
heating energy into that of gravity waves (Hines, 1965; Testud, 1970). 
Indeed, Chiu (1976) has discussed the problem of thermospheric 
response to auroral heating and points out the interdependence of 
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auroral activity and the generation of an entire spectrum of wave 
perturbations, ranging from Rossby to acoustic - gravity modes. 

An additional energy input to the nighttime thermosphere is 
thought to result from the influence of thermal tides, originating in 
the lower atmosphere. From an analysis of incoherent scatter 
measurements at Millstone Hill and St. Santin, Salah and Wand 
(1974) and Salah et al. (1975) have observed a predominantly semi- 
diurnal temperature oscillation at heights near 120 km, with a 
downward phase progression. This implies an upward flow of energy 
and the dissipation of this energy at higher altitudes. An earlier 
estimate of thermospheric heating from the semidiurnal- tide, 

performed by Lindzen and Blake (1970), led to a mean upward flux 

-4 -2 

at 120 km between 3 and 4X10 W m from this source. This 
corresponds to a rate of 7.4 to 9.8X10^^ W for the total Jaeating of 
the nighttime thermosphere from the semidiurnal tide. 

The values for the energy deficit deduced from Table 5 are 
of the same order of magnitude as the above estimates and indicate 
a need for an additional energy input to the low and mid latitude 
nighttime thermosphere. 

5. 3 Recommendations for Further Experimental Work 

It is clear that further studies of this type will require an 
improved and expanded body of neutral density and temperature 
measurements at altitudes near 120 km. Regarding the incoherent 
scatter technique, measurements made over a 24 hour period would 
be an enormous advantage even if a complete sequence cannot be 
obtained. A number of observations over several nights during 
times when the electron density is sufficiently large would be 
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valuable in improving estimates of the odd harmonics in the density 
and temperature fluctuations. Secondly, since meaningful results 
can only be obtained on a statistical basis, a regular and systematic 
program of density and temperature measurements is required to 
ascertain seasonal variations and reduce the statistical Tincertainties. 
A larger body of temperature measurements would also be useful 
in refining the conduction estimates. 

Although the emphasis here has been placed on improved 
boundary values, there remains a need for additional morning- 
evening rocket vapor trail measurements in the lower thermosphere. 
Additional twilight flights at a variety of latitudes and for different 
seasons would be of great help in revealing the nature of the wind 
system in the lower thermosphere and improving the overall heat 
balance. 

Finally, a coherent series of high latitude rocket measure- 
ments, organized in terms of geomagnetic coordinates, would be 
invaluable in the study of the wind system at high geomagnetic 


latitudes. 
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APPENDIX A: BOUNDARY CONDITIONS FROM ARECIBO 

AND ST. SANT IN 

Tables A.l and A, 2 list the temperature parameters and 120 
km densities used in this study. The results in Table A.l were 
derived from 18 days of observations at Arecibo between August, 
1974 and November, 1976. The data in Table A. 2 were obtained 
from the St. Santin facility between August, 1969 and June, 1970. 



Table A.l. Arecibo 



Aug. 12, 1974 Summer 6.40 1.8 97 264 


6.77 266 

7.14 300 

7.57 315 

17.09 1.0 304 

17.46 307 

17.83 310 

18.26 284 


779 

772 

781 

785 

922 

967 

986 

920 


Aug, 14, 1974 Summer 5.85 1.7 91 357 


6.28 334 

6.65 335 

7.02 319 

17.11 0.7 367 

17,48 366 

17.85 354 


795 

758 

850 

809 

875 

909 

972 


Sept. 17, 1974 Summer 6.29 2,0 99 364 


6,65 387 

7.09 359 

7.46 337 

7.83 339 

15.75 0.7 347 

16.12 330 

16.55 309 

16.92 277 

17.29 291 

17.72 297 


706 

723 

728 

727 

760 

964 

996 

1004 

998 

979 

1012 


t Average of three values preceding measurement. 


Conditions 


120 km Densities (M 


n(02) 


n(N2) 


.02188 

5.51x10^® 

3.86x10^^ 

1.72x10 

.04685 

4.80 

3.36 

9.92x10 

.04424 

4.68 

3.27 

8.49 

.06208 

3.56 

2.49 

7.06 

.02338 

4.37 

3.06 

1.17x10 

.02950 

6.68 

4.67 

7.99x10 

.01774 

3.12 

1.26x10 ' 

2.19 

1.12x10 

.02678 

8.81 

1.04 

.01789 

1.17x10^^ 

8.21 

9 . 15x10 

.02550 

2.65 

1.86 

9.62 

.01234 

2.11 

1.47 

1.35x10 

.02497 

4,08 

2.85 

9.15x10 

.03753 

2.54 

1.78 

6.79 

.02803 

5.83 

4.08 

6.89 

.01412 

5,85 

4.09 

1.03x10 

.04931 

8.02 

5.62 

1 . 00x10 

.02182 

6.82 

4.78 

1.26 

.04618 

7.00 

4.90 

9.97x10 

.04892 

7,37 

5.16 

1.12x10 

.03832 

6.30 

4.41 

1.09 

.03563 

5.38 

3.76 

7.53x10 

.02502 

5.30 

3.71 

9.17 

.02114 

7.13 

4.99 

1,13x10 

.02326 

8.53 

5.97 

1.25 

.02593 

5.89 

4.12 

1.06 

.01578 

8.52 

5.97 

1.46 


17 

16 


17 

16 

17 


16 

17 

16 


17 

17 

16 

17 

16 

17 
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Table A.l. (Continued) 



Sept. 18, 1974 Fall 


April 15, 1975 Spring 


6.31 

0.2 

99 

379 

755 

.04023 

7.15x10^^ 

5.00x10^^ 

6.28x10^^ 

6.68 



362 

775 

.03182 

4.61 

3.23 

7.13 

7.05 



349 

755 

.03873 

5.80 

4.06 

7.82 

7.48 



317 

772 

.04124 

5.59 

3.91 

8.58 

1.04x10::' 
8.55x10 ° 

7.85 



322 

874 

.01843 

4.65 

3.26 

8.23 



295 

834 

.03720 

6,64 

4.65 

14.91 

3.7 


444 

1020 

.01677 

2,94 

2.06 

8.91 

15.32 



424 

989 

.03374 

3.84 

2.72x107° 
4.25x107^ 
3.59x10 ^ 

2.69 

1.90x10::' 
2.98xl0r' 
2.51x10 ' 

6.83 

1.18x107' 

8.57xl0-*-° 

15.82 



421 

1208 

.00872 

16.58 



410 

1114 

.01527 

16.95 



439 

1121 

.01485 

7.32 

6.21 

2.2 

69 

407 

879 

.00922 

1.62x1oH 

7.34x10::° 

9.28x107° 

5.33x10 

T ft 

1.13x10::° 

5.i4xio:::J 

6.50x107^ 

3.73x10-*-' 

9.55x10^® 

6.62 



386 

740 

.02383 

8.75 

1.07x10 

7.02 



370 

831 

.01367 

7.47 



304 

714 

.02860 

1.54 

7.88 



321 

753 

.01705 

5.30 

3.71 

1.85 

8.32 



313 

754 

.01958 

5.12 

3.58 

1.81 , 
8.00x10 

15.43 

1.9 


378 

845 

.04334 

4.43 

3.10 

15.83 



416 

849 

.02287 

2.94 

2.06 

9.11 
1. 23x10 

16.29 



430 

1233 

.00522 

3.79 

2.65 

16.69 



412 

910 

.01219 

5.08 

3.55 

1.16 

17.09 



401 

1315 

.00472 

6.60 

4.62 

1.42 

9.54x10 ° 

17.55 



387 

922 

.01589 

7.27 

5.09 


t Average of three values preceding measurement 
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Table A.l. (Continued) 


Date 


April 16, 1975 


May 22, 1975 


Season 


Spring 


Siioraier 


Temperature Parameters 


L.S.T. 


F 

10.7 

T 

' 120 

T 

00 

S 

6.15 

2.2 

69 

432 

769 

.00819 

6.56 



396 

685 

.02089 

6.96 



389 

902 

.00826 

7.41 



406 

863 

.00797 

7.82 



397 

778 

.01846 

8.22 



366 

755 

.03768 

15.45 

1.0 


347 

932 

.01457 

15.85 



353 

891 

.02396 

16.31 



368 

961 

.01244 

16.72 



359 

882 

.01882 

17.11 



370 

969 

.01165 

6.21 

2.2 

69 

385 

739 

.02725 

6.61 



393 

914 

.00792 

7.01 



369 

787 

.02852 

7.60 



363 

788 

.01902 

8.00 



360 

794 

.03745 

8.40 



352 

810 

.02025 

15.30 

3.0 


347 

910 

.03297 

15.70 



350 

880 

.03107 

16.56 



408 

869 

.02728 

16.96 



409 

956 

.01192 

17.42 



402 

953 

.01273 


120 km 

Densities 

(M"^) 

nCOj) 

n(N2) 

n(0) 

1. 18x10^ J 
3.02x10-^^ 

8.23x10^^ 

1.25x10^^ 

2.11 

1.18 

5.78 

4.05 

1.21 

3.89 

2.73 

1.29 _ 

9.32x10'"'^ 

5.35 

3.75 

6.65 

4.66 

8.30 - , 
1.37xl0r^ 
9.45xl0r, 

1.22x10 

4.97 

3.48 

5.07 

3.55 

4.59 

3.21 

5.23 

3.66 

1.08 

3.77 

2.64 

1.18 

4.32 

3.02 

7.05x10^5 

1.04x10::^ 

6.42x10 

4.25 

2.98.. 

5.52 

3.86 

3.78 

2.64 

9.33 

3.70 

2.59 

6.28 

3.76 

2.63 

9.61 

3.93 

2.75 

8.38 

3.16 

2.21 

9.31 

4,92 

3.45 

7.65 

1.05x10-^^ 

3.59 

2.51 

4.03 

2.82 

1.00 


t Average of three values preceding measurement 
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Table A.l. (Continued) 


Date 


June 10, 1975 


June 11, 1975 


Season 


Summer 


Summer 


_3 

Temperature Parameters 120 km Densities (M ) 


L.S.T. 


^10.7 

. T 

120 

T 

CO 

S 

n(02> 


n(0) 

5.58 

0.8 

67 

350 

615 

.04911 

1. 32x10^ J 
7. 12x10 

9.27x10^^ 

9.54x10^^ 

5.98 



368 

629 

.03876 

4.98 

8.89 

l.O4xl0|^ 

8.26x10-^° 

6.38 



356 

715 

.01592 

3.68 

2.57 

6.84 



361 

746 

.01765 

3.30 

2.31 

7.24 



355 

771 

.01443 

2.30 

^ • 6 X 

9.69 

7.64 



324 

776 

.02247 

4.38 

3.07 

8.19 

16.11 

0.9 


402 

902 

.01595 

5.00 

3.50 

7.89 

16.61 



344 

947 

.01372 

5.05 

3.54 

1.09x10::' 

6.39x10 

17.01 



329 

882 

.03607 

5.20 

3.64 

17.41 



330 

806 

.03823 

5.32 

3.73 

1.97 

18.36 



334 

858 

.02154 

6.25 

4.38 

8.40 

5.73 

1.7 

66 

404 

710 

.01685 

7.57x10^^ 

5.30 

8.18x10^^ 

6.13 



393 

725 

.01755 

5.22 

3.66 

8.12 

7.00 



393 

755 

.02495 

7.47 

5.23 

6.02 

7.40 



402 

738 

.03231 

4.62 

3.24 

5.71 

7.80 



411 

866 

.01294' 

4.84 

3,39 

6.32 

16.23 

3.2 


383 

837 

.03414 

4,03 

2.82 

8.25 

16.63 



371 

935 

.01696 

4.05 

2.83 

9.88 

17.03 



350 

911 

.03108 

3.78 

2.65 

7 35 

1.03x10^ 

17.49 



356 

857 

.02224 

5.88 

4.12 

17.89 



328 

828 

,02164 

4,85 

3,39 

1.05 

18.29 



335 

943 

,01749 

5,96 

4.18 

1.01 


t Average of three values preceding measurement 
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Table A.l. (Continued) 


Date 

Oct. 14, 1975 


Jan. 20, 1976 


Season 


Fall 


Winter 


S *T « 


Temperature 

Parameters 

120 km 

Densities 

(m"^) 

F T 

10.7 120 

T 

00 

S 

n(O^) 

n(N2) 

n(0) 

7.07 

1.2 

80 414 

695 

.04568 

4.65x10^^ 

3.26x10^^ 

8.06x10^^ 

7.48 


405 

703 

.04098 

6.24 

4.37 

8.61 

7.89 


390 

719 

.04373 

9.25 

6.48 

8.59 

8.39 


383 

749 

.03292 

6.11 

4.28 


9.04 


379 

778 

.02637 

3.81 

2.67 

1.00X10-, 

9.45 


360 

799 

.03189 

7.75 7 

5.42 

9.49x10-,° 

15.04 

2.7 

348 

911 

.02989 

1. 01x10 

7.05 

1.08x10 

15.73 


349 

887 

.02651 

7.97x10^ 

5.58 

1.23 

16.13 


317 

927 

.02262 

4.87 

3.41 

1.42 

16.54 


318 

865 

.02837 

9.19 

6.43 

1.43 

17.04 


303 

943 

.01648 

7.08 

4.96 

1.91 

6.64 

2.6 

77 399 

705 

.01591 

1. 15x10 

8.04 

1.20x10^ 

7.90 


363 

698 

- .01906 

4.10x10 

2.87 

1.55 

8.30 


352 

658 

.03828 

5.48 

3.83 


15.11 

1.9 

321 

846 

.03970. 

3.53 

2.47 

9 .63x10 

15.51 


342 

829 

.05690 

3.02 

2.11 

7.89 

15.91 


387 

813 

.04253 

3.34 

2.33 

7 .69 

16.37 


384 

810 

.04884 

4.06 

2.84 

7.30 

16.77 


410 

837 

.02063 

5 . 44 

3.81 

8.74 

17.17 


441 

956 

.01021 

7.02 

4.91 

8.64 


t Average of three values preceding measurement. 
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Table A.l. (Continued) 


Temperature Parameters 


Date 


Season L.S.T. 


Jan. 21, 1976 Winter 


Mar. 17, 1976 Spring 


n(Oj) 


120 km Densities (M 
o) nCN^) n(0) 


6.81 

3.5 

75 

451 

765 

.00964 

5.42x10^^ 

3.79x10^^ 

1.07x10 

8.08 



397 

718 

.02208 

3.93 

2.75 

1.17 

8.48 



381 

725 

.02844 

5.34 

3.73 

1.10 

14.57 

3.7 


332 

839 

.06093 

4.72 

3.31 

9.21x10 

14.97 



416 

876 

.02910 

3.14 

2.20 

7.83 

15.43 



413 

875 

.02446 

3.46 

2.42 

8.65 

15.83 



411 

999 

.01138 

3.50 

2.45 

1.09x10- 

16.23 



399 

902 

.01957 

3.35 

2.34 

9.42x10 

6.32 

3.1 

75 

432 

686 

.04829 

7.11 

4.98 

8.00 

6.72 



453 

711 

,02206 

8.62 . 

6.04 

8.85 

7.18 



431 

762 

.01438 

1.02x10::' 

7.15 

1.06x10 

7.58 



416 

790 

.01474 

7.24xlo|; 

5.07 

1.07 

7.98 



385 

786 

.02764 

1. 00x10 

7.03 

8.70x10 

14.84 

2.8 


393 

930 

.01745 

5. 00x10 

3.50 

1 . 10x10 

15.25 



392 

872 

,02752 

4,62 

3.24 

9.62x10 

15.65 



396 

934 

.02122 

5.27 

3.69 

9.00 

16.10 



356 

892 

.03208 

6.43 

4.50 

9.28 

16.51 



362 

901 

.02821 

7.36 

5.15 

9.22 

16.91 



339 

851 

.02866 

5.41 

3.79 

1.18x10 


t Average of three values preceding measurement. 
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Table A.l. (Continued) 



May 13, 1976 Summer 


Sept. 22, 1976 Fall 


5.80 

1.6 

72 

408 

773 

.01250 

7,92x10^^ 

5.54x10^^ 

8.99x10 

6.21 



365 

772 

.02893 

5.35 

3.75 

6.41 

6 .61 



369 

823 

.01888 

7.03 

4.92 

7.07 

7.07 



353 

805 

.02810 

5,97 

4.18 

6.51 

7.47 



321 

758 

.03368 

8.26 

5.78 

8.88 

7.87 



295 

766 

.04037 

7,85 

5.50 

9.29 

15.15 

1.9 


1 420 

904 

.01900 

5.00 

3.50 

8.43 

15.55 



1 395 

974 

.01436 

5.42 

3.80 

9.62 

16.41 



372 

909 

.01459 

5.37 

3.76 

1.26x10 

16.81 



352 

893 

.01928 

5.52 

3.86 

1.12 

17.27 



346 

876 

.01997 

6.57 

4.60 

1.15 

6.63 

3.5 

70 

323 

816 

.01341 

5.78 

4.05 

1.80x10 

7.02 



341 

747 

.02066 

3.90 

2.73 

1.45 

7.49 



335 

717 

.02817 

4.06 

2.84 

1.46 

7.89 



320 

725 

.04823 

5,09 

3.57 

1.18 

8.28 



335 

746 

.04226 

3.93 

2.75 

1.09 

8.74 



326 

746 

,04444 

3.63 

2.54 

1.18 

15.26 

3.2 


432 

780 

.04320 

2.62 

1.83 

9.47x10 

15.66 



441 

839 

.02716 

2.60 

1.82 

8.58 

16.05 



436 

858 

.02270 

3.60 

2.52 

8,84 

16.52 



428 

885 

.01892 

5.41 

2.79 

9.16 

16.92 



401 

897 

,02466 

4,94 

3.46 

8.13 


17 

17 

16 


N5 


t Average of three values preceding measurement . 



Table A.l. 



Oct. 19, 1976 Fall 


6.41 

0.8 77 

384 

717 

6.81 


376 

687 

7.68 


362 

693 

8.08 


372 

697 

8.48 


367 

748 

15.64 

2.8 

376 

877 

16.04 


366 

868 

16.90 


372 

897 

17.30 


356 

864 


Nov. 16, 1976 Winter 7.52 1.0 73 341 


8.10 328 

8.58 356 

9.03 365 

9.48 338 

14.98 1.7 347 

15.38 346 

15.91 328 

16.49 352 

16.89 336 

17.29 321 


672 

708 

668 

767 

735 

899 

950 

892 

856 

879 

861 


t Average of three values preceding measurement 


inued) 




120 km Densities (M 


nCN^) 


.01899 6.54x10^^ 

.03214 5.01 

.03832 3.72 

.03338 2.91 

.02298 3.57 

.02993 3.86 

.03064 5.30 

.01778 4,99 

.02289 3.48 

.05455 6.25 

.05666 3.26 

.06963 3.13 

.01954 2.98 

.04490 3.43 

.02140 3.91 

.02081 3.69 

.02692 5.61 

.02862 3,43 

.02408 4.28 

.02483 6,47 


4.58x10^^ 

1.08x10^ 

3.51 

9.97x10- 

2.60 

1.04x10' 

2.04 

1.10 

2.50 

1.13 

2.70 

1.01 

3.71 

1.06 

3.49 

1.30 

2.44 

1.26 

4.38 

1.05x10 

2.29 

9.48x10 

2.19 

8.69 

2.08 

1.20x10- 

2.40 

1.03 

2.73 

1.18 

2.58 

1.02 

3.93 

1.08 

2.40 

1.00 

3.00 

1.10 

4.53 

1.21 


17 

16 

17 


17 

16 


17 
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Table A.l. (Continued) 


_3 

T emperature Parameters 120 km Densities (M ) 

Date Season L.S.T. Kp^ ^10.7 ^120 ^ 


Nov. 17, 1976 Winter 6.90 

2.2 

75 392 

753 

.01797 

4,17x10^® 

2.92x10^^ 

1.07xloJ-J 

7,37 


360 

711 

.06276 

3.18 

2.23 

8.08x10-^'’ 

' 7.77 


370 

736 

.05167 

3.55 

2.49 

7.98 

8,16 


369 

721 

.04695 

3.90 

2.73 

9.37 

9.01 


396 

771 

.02755 

3.39 

2,37 

9.22 

9.80 


377 

777 

.02312 

4.35 

3.05 

1.21x10 

14.96 

1.3 

336 

866 

.02804 

3.74 

2.62 

1.11 

15.36 


330 

858 

,02640 

4.65 

3.26 

1.21 

15.83 


314 

844 

,03667 

4.32 

3.02 

1.10 

16.23 


318 

821 

.04185 

5,05 

3.54 

1.07 

16.62 


364 

1000 

,00973 

2.32 

1.03 

1.54 

17,08 


311 

929 

,01808 

3.27 

2.29 

1.33 


t Average of three values preceding measurement 
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Table A. 2. St. Santin Boundary Conditxons 








Temperature 







L.S.T. 

SF/SS 

K-' 

P 


Parameters 

120 

km Densities 

(M 

Date 

Season 

F 

10.7 

T 

120 

T 

CO 

n(02) 

nCN^) 

n(0) 

Aug. 

26, 1969 

Fall & 

5.30 

1.2 

153 

398 

1025 

4,llxlo}'^ 

7.34x10 

2.88x10?-7 

5.13x10 

9.26x10^^ 

8.42x10 



Summer 

18.70 

2.4 


428 

1170 

Oct , 

7, 1969 

Fall 

6.36 

2.5 

145 

460 

820 

4.86x10^^ 

4.48x10 

3.40x10?-^ 

3.14x10 

1.41x10^^ 

7.76x10 




17.64 

0.6 


485 

1115 

Nov. 

25, 1969 

Winter 

7.46 

1.1 

191 

455 

895 

4.23x10^^ 

7.65x10-^^ 

2.96x10^7 

5.35x10 

1.08xloj-^ 

1.57x10-^^ 




16.54 

1.8 


357 

1125 

Dec . 

23, 1969 

Winter 

7.69 

2.0 

162 

422 

860 

4.91x10?-^ 

9.87x10 

3. 44x10 ?- 7 
6.91xl0-^ 

1.59x10^^ 

1.24x10-^' 




16.31 

2.9 


380 

1100 

Feb. 

5, 1970 

Winter 

7.09 

2.4 

125 

300 

850 

2. 95x10 
4.19x10 

2.07x10^7 

2.93x10 

2. 87x10 
1.17x10 




16.91 

2.1 


420 

1050 

Mar . 

17, 1970 

Spring 

6.09 

0.6 

135 

398 

960 

6.87xloj-^ 

9.71x10 

4 .8lxlo|? 
6.80x10 

17 

1.14x10, ^ 
7.52x10 




17.91 

2.1 


490 

1150 

Apr. 

28, 1970 

Summer 

5,04 

2.1 

147 

408 

970 

6. 21x10 
4.75x10 

4.35x10^7 

3.33x10 

1.26x10^^ 

9.27x10 




18.96 

1.0 


475 

1105 

June 

9, 1970 

Summer 

4 . 36 

1.2 

134 

443 

1060 

6.95x10^5 

6.78x10 

4.87x10?-^ 

4.74x10-^^ 

8.35x10?-^ 

7.20x10 




19.64 

1.6 


457 

1170 

Sept 

. 1, 1970 

Fall 

5.45 

1.9 

151 

425 

1160 

6.05x10^^ 

7.02x10 

4.23x10?-^ 

4,92x10 

7.86x10^^ 

8.86x10 




18.55 

4.0 


408 

1185 


+ 


Average of three values preceding measurement. 
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fhoro haa 2>con conalietabXe spocuJAtion as to the onorgy balaaco of cAo 
thoroospAoro And in particuiar aAoue tAo fraction of cAo cotaJ energy input auppliod 
by uitra^vioioe radiation The problon is conaidoraMy aimplified Ay eonsidoriny tAo 
energy halanco of tAo niyAttimo Ao'niapAoro alone Sunrise and sonaet vapor trail 
n>ea%arcinonts provide data on tAo wind systena at tAo corninator boundary^ and 
torporaturo meaauremonta provide inforaiation on the vertical oftoryy conduction 
North'‘South hinds fron high latitude vapor traii roaaurononca provide a loeasaro of 
the energy input from auroral ptoeossos The discropanctj between the inputa froai 
tAo day Aoaiiop^ioro and the observed loss rates is diocuaood in tornu of ion«noutrai 
processes and gravity wave inputs 
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TAoro has been conaidoraAio speculation as to tho energy balance of tAo 
tAormoapAoro and in particoiar about tAo fraction of tho total onoryy input supplied 
by ultra-violet radiation TAo preblen is considerably sinpjifiod by considering tAo 
energy bilanco of tAo niyAttimo hemisphere aiono 5unriao and sunaot vapor trail 
rnoaauronones provide data on cAo »/ij)d oyotona at tAo torminator boundary, and 
tonporacurc n»aooromoncs provide inforr*ation on tAo vorticai onoryy conduction 
//ortA'OOutA vinds from AiyA iatieudo vapor traij tnaasureaenta provide a toeasure of 
cAo energy input from auroral processes The discrepancy between the inputs from 
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tAormospijoro and in particular about the fraction of tho total energy input supplied 
by uitra-vJo2ofr radiation The jarobJem is considerably sitopllficd by considering tho 
energy balance of the nighttino Aontisphoro aionc Sunrise and sunset vapor trail 
laoasurencnts provide data on tAo wind syatoms at tAo torminator boundary, and 
temperature noasuronents provide infornation on tho vortical onorgy conduction 
NOrth’SOUth winds from high latitude vapor trail Bteasureaants provide a moaeuro of 
eJ>o onorgy input froo auroral processes The discropanty between the Inputs from 
the day hemisphere and the observed loss rates is discussed in terms of ion-neutral 
processes and gravity wave inputs 
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There Aas boon considerable speculation as to tho energy baianco of tAo 
cAermospAoro and in particular about the fraction of tho total energy input suppiied 
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energy balance of tho nighttime homlaphore alone Sunrise and sunsat vapor trail 
moasuromonts provide dati on tho wind aystens at tho terminator boundary, and 
tomporatuzQ racasuramonts provide infornation on tAo vertical energy conduction 
PortA-soucA winds from high latitude vapor trail moasuromonts provide a oeasuro of 
tho energy input from aurorai processes TAo discrepancy between the inputs from 
the day Aonisphero and the observed loss rates is discussed in terms of ion*neutrai 
processes and gravity wave inputs 
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by ultra-violet radiation The probJen is ccnsidorablg aioipJifJcd by conoidering the i p Etructuro of the Upper 
energy balance of tho nigbttJrio honisphore alone Sunrise and sunset vapor trail Ataosphoro 

noasuranonta provide data on the wind ayatens at tho tominacor boundary, and 19 3 Uetoroapbore 

tomporaturo neasuxomonta provide information on the vertical energy conduction 194 Interaction Datveon tho 

North-south winda from high latitude vapor trail moasuromonts provido a moosaro of Neutral Atoospitore and 

the energy inpot from auroral procossos The discrepancy between the inputs from tiie Jonoaphoro 

tho day homlsphoro and tho observed loss rates ia discussed in terms of ion-neutral 
processes end gravity wave inputs 
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There has boon conaiderablo speculation as to tho energy balance of the 
thermosphere and in particular about the fraction of the total energy input supplied 
by v,2bra-iJeJet radiation Tho probJon is considerably siaolifiod by considering the 
energy balance of tho nighttime hemisphere alone Sunrise and sunset vapor trail 
measurements provide data on the wind systems at tho terriinator boundary/ and 
tertporacuro moasuremonts provide information on the vertical cnorgy conduction 
Worth-south winds from high JatJtudo vapor trail moaeuroments provide a auasaro of 
the cnorgy input from auroral procoaaos Tho discrepancy between the anputs from 
the day homisphero and the observed loss rates ia discussed ia terms of ion-neutral 
processes and gravity wave inputs 
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Glonar, David Alan, The Energy Dalanco of the Pighetimo Thermosphere, The psv-iRL-sci-dS? 

Ionosphere Research Laboratory, Electrical Engineering East, ' 

University Pirh, Pennsylvania, 16802, 1977 classification Numbers 

There has boon Considerable speculation as to the energy balance of the 15 2 E-Region 

thcrrtosphQzo and in particular about tho fraction of the total energy input supplied 1 5 J F-Rogion 

by ultra-violoc radiation Tho problom ia conaidorablg simplified by considering the 1 9 Structure of tho Dppor 
energy balance of tlte nighttimo homlaphoro alone Sunrise and sunset vapor trail Atmesphoro 

masurooants provido data on tha wind systems at tho terminator boundary and 19 3 Uotoroaphoro 

temperature measurements provido Inforaation on the vortical energy conduction 124 interaction between tho 

Porth-soutii winda from high latitude vapor trail moaaurcmonta provide a neasure of Neutral Attoosphoxo and 

tha energy input from auroral procossos Ti>e rfiscroponcy botwoon the inputs from tho lonoapltora 

tho day hemisphere and tho observed less rates is discussed in terms of ion-neutral 
processes and gravity wave inputs 
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